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PREFACE 


The  1995  Articulated  Total  Body  (ATB)  Model  Users’  Colloquium  was  held  at  the  Wright- 
Patterson  Air  Force  Base,  Dayton,  Ohio  on  June  13  and  14, 1995.  This  Colloquium,  sponsored  by 
the  Armstrong  Laboratory  (AL),  U.  S.  Department  of  the  Air  Force,  brought  together  more  than 
seventy  users  of  the  ATB  model  and  its  derivatives  (CVS,  Cal- 3D,  and  DYNAMAN).  The  two-day 
conference  offered  the  opportunity  to  present  and  exchange  the  latest  ATB  modeling  techniques  and 
applications.  Dr.  Thomas  J.  Moore,  Chief  of  the  Biodynamics  and  Biocommunications  Division  of 
AL  welcomed  the  attendees.  Eighteen  invited  presentations,  group  discussions,  and  interactive 
exercises  covered  areas  like  model  algorithms,  harness  belt  and  airbag  modeling,  data  base 
development,  vehicle  and  aircraft  crashworthiness,  dummy  and  real  human  modeling,  and  accident 
reconstruction  applications.  A  formal  ATB  Users’  Group  was  created  to  organize  regular  meetings, 
to  set  up  bulletin  board  for  exchanging  information,  and  to  offer  other  opportunities. 

As  this  colloquium  was  designed  to  be  an  informal  meeting,  please  do  not  reference  the 
presentations  appearing  in  these  proceedings  in  the  open  literature. 


The  Organizing  Committee  would  like  to  thank  AL  and  Systems  Research  Laboratories,  Inc.  for 
their  support  of  this  colloquium.  Many  thanks  are  extended  to  the  speakers  for  investing  the 
resources  needed  to  prepare  the  presentations  which  are  contained  in  this  proceeding. 
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SYSTEM  OF  CONNECTED  Rl&lD  BODIES 


External  forces 


FREE  BODY  CONFIGURATION 


FORCES 
TORQUES 
LINEAR  JOINT 
ANGULAR  JOINT 
SPECIAL  CONSTRAINTS 
FLEXIBLE  ELEMENT 


SYSTEM  EQUATIONS 


NOMENCLATURE 

A-,  Bjj  =  MATRICES  USED  IN  SYSTEM  EQUATIONS 
f  =  JOINT  CONSTRAINT  FORCE 
M  =  MASS  MATRIX 
q  =  SPECIAL  CONSTRAINT  FORCE 
t  =  JOINT  CONSTRAINT  TORQUE 
U1  =  EXTERNAL  FORCE 
U2  =  EXTERNAL  TORQUE 
Vj  =  CONSTRAINT  FORCING  FUNCTION 
7T  =  FLEXIBLE  ELEMENT  CONSTRAINT  TORQUE 
$  =  INERTIA  MATRIX 
X  =  TRANSLATIONAL  ACCELERATION 

* 

CO  =  ANGULAR  ACCELERATION 

Figure  2  SYSTEM  EQUATIONS 
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FIGURE  2. 
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BELT  1 


Figure  3  BELT  HARNESS  MODEL 


JUNCTION  POINT 


Figure  4  REFERENCE  POINT  SELECTION 
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Figure  6  POINTS  OF  A  SAMPLE  BELT  SYSTEM 
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Figure  7  MATRIX  FORM  OF  CONSTRAINT  EQUATIONS 
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Figure  1  FIFTEEN  SEGMENT  MAN-MODEL 
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STANDARD  15- SEGMENT  OCCUPANT 


Presentation  of 


David  J.  Biss 

Automotive  Safety  Analysis  Corporation 

at 

The  1995  ATB  Model  User's  Colloquium 
June  13,  1995 


USE  OF  ATB  IN  A  REAR  IMPACT  RECONSTRUCTION 


Abstract 


An  interesting  field  case  study  was  undertaken  involving  a  180  deg.  rear 
impact  in  which  the  driver's  seat  back  did  not  appreciably  deform  or  collapse 
rearward  raising  an  issue  regarding  the  efficacy  of  a  head  restraint  adjusted  in  the 
full  down  position.  A  chosen  component  of  the  case  study  was  the  use  of  the 
ATB  model  to  analyze  occupant  kinematics  and  dynamics  with  alternative  height 
head  restraints.  A  methodology  will  be  presented  to  the  ATB  User  s  Colloquium 
showing  the  development  of  the  case  occupant  anthropometry  and  the  rendering 
of  that  dimensioned  occupant  with  a  GEBOD2  generated  surrogate  for  ATB  work. 
This  ATB  study  is  still  in  progress  and  questions  will  be  presented  to,  and  answers 
solicited  from,  the  very  knowledgeable  body  of  attendees  at  the  ATB  User's 
Colloquium  to  elicit  some  insight  regarding:  1)  Assisting  the  diagnosis  of  some 
ongoing  specific  user  problems,  and  in  the  further  successful  application  of  ATB  to 
this  rear  impact  accident  mode  study;  and,  2)  Presenting  to  the  ATB  User  s 
Colloquium  a  user's  case  study  demonstrating  the  need  for  an  Internet  or  other  on¬ 
line  forum  so  that  the  considerable  collective  insights  gained  by  the  very  diverse 
group  of  users  during  their  many  varied  applications  of  the  model  would  be 
available  to  others  with  similarly  recurring  user  problems. 
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Seat  Structure  Modeling 

and 

Fixing  a  Bug  in  the  ATB  Harness-Belt 
Unloading  Algorithm 


Ed  Sieveka 

U  n  iversity  of  Vi  rg  ini  a 
Automobile  Safety  Lab orato ry 


ATB  Simulation  and  Automobile  Safety 
Research  at  the  University  of  Virginia 
Automobile  Safety  Laboratory 

■  UVA  has  been  working  with  ATB  for 
about  6  years 

■  ASL  experiments  provide  a  ready 
source  of  ATB  validation  data 

■  ATB  is  used  in  conjunction  with 
research  at  the  ASL  for  experiment 
design,  analysis  of  results,  and 
parametric  studies 
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The  UVA  Deceleration  Sled:  Front  View 
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Toepan  Intrusion  Using  the  UVA 

Compound  Pendulum 


The  UVA  Head-drop  Device 


WMiliilMl 


Proceedings  of  the  1995  ATB  Model  Users’  Colloquium 


Part  I:  An  Unloading  Problem  in 
the  ATB  Harness  Belt  Routine 

■  Experienced  excessive  occupant  rebound 
in  Simula ti on  of  a  frontal  impa ct  s I e  d  te st 

■  Reduced  harness-belt  elasticity  via  the  R 

'  and  G  factors 

■  No  significant  improvement  in  occupant 
kinematics  occurred 

■  Shoulder-belt  load  time-histories  showed 

'  absence  of  expected  inelasticity 
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Hypothetically  Equivalent  Systems  to  Test 
ATB  Integration  Step-Size  Sensitivity 
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However,  reducing  the  step-size 
in  the  occupant  simulation  had 
no  effect  on  the  results. 
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Further  examination  of  the  belt  load  time-histories 
revealed  that  inelastic  unloading  at  the  upper 
shoulder  belt  segment  did  not  be^in  until  thi  load 
fell  below  the  maximum  load  achieved  by  the 
lowest  segment. 


The  unloading  problem  is 

caused  primarily  by  a 

failure  of  ATB  to  store 
load  history  information 
for  all  but  the  last  segment 


3-8 
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The  problem  is  fixed  by  ... 

■  Increasing  the  dimensions  of  array 
NTHRNS  from  (20)  to  (20.25)  to 
provide  storage  for  each  belt  point 

■  Store  load  history  information  for 
each  belt  point  in  the  TAB  and  NTAB 

arrays 

■  Determine  unloading  for  each  belt 
point  separately 


Proceedings  of  the  1995  ATB  Model  Users'  Colloquium 


Sphere  Velocity  Time-Histories:  Original 


There  is  still  a 
noticeable  sensitivity 
to  integration 
step-size  in  the 
harness-belt 
algorithm,  but  the 
magnitude  of  the 
effect  is  con siderably 
reduced 


Sphere  Velocity  Time-Histories:  Corrected 
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Part  II:  Using  ATB  to  Model  the 
Structure  of  a  Passenger  Car  Seat 

■  Performed  in  support  of  an  NHTSA 
project  to  develop  a  MADYMO  seat 
model 

■  Useful  for  studying  seat-back 
response  in  rear  impacts  or  in  frontal 
impacts  with  integrated  seat  belts 

■  Allows  loading  information  to  be 
obtained  with  much  greater 
efficiency  than  with  a  F.E.  model 
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11  MPH 
Rear  Impact 
Seat  Test  on 
UVA  Sled 


Modeling  Considerations: 

i  Seats  are  not  symmetric  side-to-side 

nor  are  the  applied  loads 
i  Major  seat  load  paths  and  degrees  of 

freedom  must  be  simulated 

i  The  tree-structured  nature  of 

MADYMO  and  ATB  must  be  worked 

around 
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Method  I:  Double  Horseshoes 
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Spring-Damper  Torsional  Element  for 
Seat  Back  Hinge  Joint  Simulation 


Method  II:  Coincident,  Constrained 
Segments  Simulating  Closed  Loops 
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Results  of  30  mph  frontal 
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NEW  USTMINT BELT 
MOBIL 


Louise  Obergefell 
Annstrong  Laboratory 
Wright-Patterson  Air  Force  Base,  Ohio 


CURRENT  BELT  MODELS 

•  Individual  Spring  Elements 

•  Interacting  Spring  Elements 

-  Stress  in  one  spring  affects  other  springs 

•  ATB  Harness 

-  Series  of  belt  segments 

-  Endpoints  move  on  one  ellipsoid  surface 
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OBJECTIVE 

Develop,  Implement,  &  Validate 
Restraint  Belt  Model 

•  Applies  Forces  To  The  Body 

»  Uses  Measured  Material  Properties 

•  Allows  Slippage  Of  Belt  On  Body 

•  Includes  Multiple  Belts  & 

Retraction  Systems 
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POINT  CONTACT 
WITH  TWO  ELLIPSOIDS 


L\  " 

/\\ 


VAUDAXIQN 


•  Elementary  Simulations 

-  One  Free  Point  Belt 

-  Multiple  Free  Point  Belt 

-  Buckle 

-  Multiple  Contact  Ellipsoids 

•  Full  Body  Simulations 

-  Sled  Test  Simulations 
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THREE  POINT  BELT 
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THREE  POINT  BELT 


SEVEN  POINT  BELT 

OFFSET 
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Belt  Stress  Force  (lb) 


CONCLUSIONS 

*  Restraint  Belt  Model  Provides  Reasonable 
Response 

-  Belt  Strain 

-  Surface  Deformation 

-  Belt  Slippage 

-  Across  multiple  surfaces 

-  Force  Application 

» Buckles  &  Other  Mechanisms 
8  User  Defined  Belt  &  Surface  Properties 

•  Converges  Quickly  in  Most  Cases 


FUTURE  EFFORTS 


•  Full  Body  Validations 

*  Check  For  Divergence  &  Modify 

Belt  Equations 

*  Method  for  Obtaining  Initial  Conditions 

•  Additional  Belt  Mechanisms 
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Simula’s  ATB  Experience 
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Cockpit  Air  Bag  Systems  (CABS) 
Development 
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Government  Products,  inc. 


CABS  Modeling  Requirements 


5-4 


Proceedings  of  the  1995  ATB  Model  Users'  Colloquium 


Government  Products,  Inc. 


Air  Bag  Deployment  Modeling 
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MSC/DYTRAN 
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MSC/DYTRAN  -  ATB  Contact  Interface 


C 

o 

o 

TJ  O 
U)  £ 

”  ^ 
CO  £ 

«  o 

TJ  = 
0  CO 


0)  CD 
T5  H 
O  < 

£  0 

if 

o  p 


5  i2 

-  c 

Z  0 

>1 

°  <t 

■E  cc 

CO  !r 
T3  >■ 

O  Q 

a® 

■0  o 


2  ° 

~  o 

"O  0 

£  C 
O  O 
0  <■> 

E  o 
0  +- 
3  "O 
o  0 

05  o 

0  = 

c  10 
0  0 

E  0  to 
2  ^  ■£ 
CD  O  o 

=  2  E 
£  5  I 

0  CO  ^ 

z?z 
<  0  jJ 

cc  c  cc 

t 


Proceedings  of  the  1995  ATB  Model  Users'  Colloquium 


Government  Products,  Inc. 


other  rigid 
ellipsoids,  shell 


CD 

C 

CD 

n 

E 

o 

E 

■o 

c 

CD 


CO 

+-» 

c 

CD 

E 

0) 

CD 


©  ^ 

;  w  S 

-  ■-  45  c  <o 

*  1 1  E  | 

5  I  til 

)  O  CD  CD  CD 


CD 

£ 

■  HI 

©  c 

■O  o 

r-  +3 

©  .E 

O  J= 


u 

9 

E 

K 


5-8 


Proceedings  of  the  1995  ATB  Model  Users’  Colloquium 


Government  Products.  Inc. 


CABS  Simulation 
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CABS  Simulation  Model 
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Crash  Pulse  Application 
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Simulation  Results 
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Summary 
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Further  Enhancements 
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Suggested  ATB  Enhancements 
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MSC/DYTRAN  Features 
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Benefits  of  ATB  in  MSC/DYTRAN 
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Benefits  of  ATB  in  MSC/DYTRAN 

(continued) 
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Integration  of  ATB  in  MSC/DYTRAN 

Methodology 
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Integration  of  ATB  in  MSC/DYTRAN 

Methodology 
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MSC/DYTRAN  input  deck 


Positioning  of  ATB  Occupant  Model 
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Visualization  of  Coupled 
A  TB-MSC/DYTRAN  Results 


cd  5 
cc  h; 

fie 


S  2  c 

r*  P  c  q> 

^  §C(0 

(b  j)  S  §  o 


cc<£ 
o  </) 

ll 

c  c/> 

o  „s> 

<2 

CO  ^  3 

|  §£■ 

0*0  3 

£ 

Q.*§  j§ 
"5  "S  -S2 

o  s-c 

Oft  O  CD 

||.i 

^  RS  -t- 


W) 

II 

o  ^ 

o  q> 

==  CD 

cc  -C 

o  (0 


CC  °  ^ 

cc  qP  3 
CO  £  O 

^  -g 

OQ  N 

CO  >*  ig 

•S*5ft 

■2  (/)  3 

q>  ^  .jo 

CQ  £  ^ 


co  vq  o: 
0)  Oijv 
^  c  O 

CC  73  CO 
K*  Q)  CC 

^r-S-c 

Q  "cc  ■§ 

o  3  5 

CO  .JO  ^2 

is  ^  o 

i  y  Q 

b  &  7c 

^  CD  g 

o  cd  cd 

■5.-°  E 

JC  c 

54  8 

OS  u 


ft 

i  «: 

11 

%  e 

ill 


Proceedings  of  the  1995  ATB  Model  Users’  Colloquium 


1995  ATB  Users’  Colloquim 


Development  of  Enhanced  50% 
HYBRID  III  ATS  Occupant  Model 
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Validation  of  Enhanced  50%  HYBRID  III 
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LEGEND 

- CALCULATED  UPPER  TORSO  ACCE  LLERATION 

— © — EXPERIMENTAL  UPPER  TORSO  ACCE  LLERATION 


ACCELERATION  TIME  HISTORIES  FOR  LOW  SEVERITY  LEVEL 


Resultant  Acceleration  (g) 
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SEVERITY  LEVEL 


LEGEND 

- CALCULATED  LOWER  TORSO  ACCELLERATI ON 

— O —  EXPERIMENTAL  LOWER  TORSO  ACCELLERATION 
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Helicopter  Crash  Characteristics 

•  Wide  range  of  pitch  and  roll  conditions 
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ATB  Model  Applications 
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Energy-Absorbing  (EA)  Seating 
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EA  Seat  Simulation  With  ATB 
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Compare  with  Validated  Model 
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Occupant  dynamic  responses  for 
Pulse  2  (contd) 
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Time  =  100  msec 


Occupant  dynamic  responses  for 
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Crash  Pulse  4 
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Occupant  dynamic  responses  for 
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Time  =  0  msec 


Occupant  dynamic  responses  for 
Pulse  4  (contd) 
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Time  =  100  msec 


Occupant  dynamic  responses  for 
Pulse  4  (contd) 
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Time  =  140  msec 


Buttocks/seat  pan  contact  force 
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Fuel  Cell  Installed  in  Helicopter 
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ATB  Model  of  Fuel  Cell 
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Predicted  Fuel  Cell  Motions 
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Proprietary  Bug  Fixes  with  Public  Domain 
Software 
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Dynamic  Mass  Property  Analysis 

—  Combining  the  6D0F  model  with  the  ATB  model 

-  Enables  injury  projection  throughout  the  escape 
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Xcg-Xcg(0) 


Moments  of  Inertia  Variation 
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The  final  analysis  enables  a  more  complete  solution 


Introduction  of  Deformable  Segments  in  The  ATB  Model 
1-  Human  and  Manikin  Neck  Models 
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OBJECTIVE 


1-  To  accurately  model  segments  with  significant  structural 
deformation  in  addition  to  the  rigid  body  motion. 

2-  To  predict  head/neck  motion  more  accurately. 

3-  To  help  design  a  more  realistic  manikin  neck. 
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APPROACH 


*  Revise  ATB  mathematical  formulation  to  incorporate 
vibrational  deformation  modes  of  individual  segments 
in  addition  to  the  rigid  body  modes. 

*  Develop  finite  element  model(s)  of  the  deformable 
segment(s)  and  determine  the  relevant  natural  frequencies 
and  mode  shapes. 

*  Create  a  data  file  for  each  deformable  segments  containing 
nodal  positions,  natural  frequencies,  and  mode  shapes. 

*  Revise  the  ATB  “.ain”  (input)  file  to  accept  and  read  the 
deformable  segment  data. 


LIMITATIONS: 

*  Plastic  (permanent)  deformation  can  not  be  modeled. 

*  Rotational  deformation  must  be  less  than  1  rad. 
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TNPTJT  DATA 


*  There  are  two  postprocessors  written  for  ANSYS  and 
NIKE3D  to  create  the  additional  input  file(s)  needed 
for  the  ATB. 

*  Minimum  additional  data  must  be  entered  in  the  ATB 
“.ain”  file.  These  data  include: 

-  No.  of  deformable  bodies  (segments) 

-  ATB  segment  no.  of  the  deformable  body 

-  Deformable  body  data  file  name 

-  Modal  damping  ratios  (estimated) 

-  Initial  modal  positions  and  velocities  (usually  0) 

-  Joint  attachment  node  numbers 


NOTE:  The  revised  ATB  computer  program  is  compatible 
with  the  previous  ATB  programs  for  rigid  body  modeling. 
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Hybrid  III  finite  element  mesh 
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Hybrid  ni  Head/Neck  Pendulum  Test  Model 
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Time  (sec) 

Hybrid  III  rotations  in  HNP  20°  flexion  test 
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Rotation  (deg) 


Time  (sec) 

Hybrid  III  rotations  in  HNP  60°  flexion  test 
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Hybrid  III  rotations  in  HNP  120°  flexion  test 
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Rotation  (deg) 


Time  (sec) 


Hybrid  III  rotations  in  HNP  65°  lateral  test 
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Moment  (in-lb) 


0.00  0.05  0.10  0.15  0.20  0.25  0.30 

Time  (sec) 

,  Moments  for  60°  flexion  and  65°  lateral  tests 
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Torso-Head-Neck  Model  for  Sled  Test 
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Clinical  Motivation  for  Knee  Testing 


the  knee  joint 


ORL  Diarthrodial  Joint 
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Replace  SPG  data  used  in  models  with  MRI  and  CAT 
scan  data  from  clinical  studies 


Purpose/Applications 


Q-  0 
0  "O 


0  < 
“O  _ 


0 

o  £ 

0  .E 

D  O 

E  T 


•>,  .E 

1—  (— 

0 

O  0 

0  E 

o  0 
0  & 

0  Q- 

0  O 

£  0 
^  -I— ' 

>  £ 

“O  O 


T3  ^ 

0  c 
o  ■= 

o  o 

- * 

CL  0 

0  ■*-' 
o  c 

O)  ~ 

a  r 

i—  *23 

3  '  > 
0  5 

0  0 
13  C 

o  J“ 
~  0 

z:  ±15 
0  -0 

>  a 

°  o 

0  0 
Hh-»  ^3 

J>  § 

®  “ 

5  § 

5>i 

3  2 


C 

o  0  ;o 

&  s.  ® 

CD  E 

£  go 

°  2  8> 
§  »:2 
£  0  -c 

0  ”5  0 
-£  O  O 

0  E  0 

o  ^  ~ 
T-,  LU  C 
?  u.  r 

«  -o  ■= 

r-  CD  > 

.2  «  .E 

■5(5  2 

o  -o  % 

°-  ©  T3 

°r§ 

IS  8 

J  _C  J 
®  c§  « 

*■’  ®  5) 

®  2  ® 

_c  0  s- 

-  *t  © 

0)  =3  £ 

w  w  ^ 

3 


Proceedings  of  the  1995  ATB  Model  Users'  Colloquium 


ORL  Approach  to  Joint  Modeling 

Experimental 
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and  proteoglycans 

Determined  at  same  site  as  biomechanical  properties 


ORL  Approach  to  Joint  Modeling 

Experimental 


r 


0 

CL  C 

0  ■£ 

"C!  05 

0 

H— '  L- 
0  0 

C  — 
CO  >, 

S 


Q)  g 

0  t 
0  Z5 
‘O  0 

2  CD 
CL 

c  J 

■3  £ 

s  2 

o  j_ 

o 

0  .2 
CL  t 
0  CO 


©  cl 

3  | 

1=  E 

=  CO 
CO  </} 

0  © 
O  5 

.o  .2 

CD 

■E  0 

>>  CD 
2—  0 
0  r= 

-a  tr 

C  0 
13  O 


Proceedings  of  the  1995  ATB  Mode!  tiers'  Colloquium 


experimental  conditions  using  ORL  joint  testing  rig 


Spotlight 


c 


o  CO 
O-43-O 

o  (0  o 

hr  '  r- 

CL  d  h— * 

"TT  0  O 

c  E  £ 
0  ” 
^.0  0 

Ilf 

§  «■§ 

cc  5  5 

^  (O 

o  $  2 

O  >  CO 
CO  ®  c  o. 

5^,'c  £ 
t  O)  o  o 

_Q  I  ZO 


9-10 


Proceedings  of  the  1995  ATB  Model  Users’  Colloquium 


Surface  curvatures  for  stress  analysis 
Contact  areas  in  conjunction  with  kinematics 


ORL  Approach  to  Joint  Modeling 
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Kinematic  Studies 
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-  obtain  joint  kinematics 

-  obtain  related  joint  kinematics 


Data  Obtained  With 
Coordinate  Measuring  Machine 
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Data  Obtained  With 
rdinate  Measurina  Machine 
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Relationship  of  precise  SPG  cartilage  and  bone 
surfaces  to  local  coordinates 
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Relationship  of  precise  SPG  cartilage  and  bone 
surfaces  to  local  coordinates 


Data  Obtained  With 


9-18 


Proceedings  of  the  1995  ATB  Model  Users'  Colloquium 


Data  Obtained  With 
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Transform  both  sets  of  data  to  local  coordinates 
of  bone 
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Produces  more  consistent  kinematic  results 
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ATB  Model  Modifications 
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-  SPG  surfaces 


ATB-SPG  Prescribed  Motion 
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6  DOF  spline-fit  C  cards  (option  4)  used  for  displacement 
and  orientation  to  interpolate  between  known  points 
for  a  continuous  motion  with  constant  velocity 
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movie  format  using  in-house  View_SGI 


ATB  -  SPG  Surface 
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Applied  symmetric  spring  dampers  with  various 
damping  coefficents  to  observe  equilibrium  results 


SPG  Plane  Displacement 
Completely  Unconstrained 
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Time  (msec) 


Future  Directions 
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study  (n=6)  with  known  loads  and  kinematics  tor 
0°  through  120°  of  flexion 
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Relate  the  joint  coordinate  systems  to  the  principal 
axes  coordinate  system  of  the  segment 
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"Biphasic  indentation  of  articular  cartilage:  I.  theoretical  analysis,"  Mak  A.F.,  Lai  W.M.,  Mow  V.C.,  /.  Biomechanics, 

20:703-714,  1987. 

"A  finite  element  formulation  of  the  nonlinear  biphasic  model  for  articular  cartilage  and  hydrated  soft  tissue  including 
strain-dependent  permeability,"  Spilker  R.L.,  Suh  J.-K.,  Mow  V.C.,  in  Computational  Methods  in  Bioengineering, 
ASME,  ed.  R.L.  Spilker  and  B.R.  Simon,  New  York,  New  York,  pp.  81-92, 1988. 

"Biphasic  indentation  of  articular  cartilage:  II.  a  numerical  algorithm  and  an  experimental  study.  Mow  V.C.,  Gibbs 
M.C.  Lai  W.M.,  Zhu  W.B.,  Athanasiou  K.A.,  /.  Biomechanics,  22:853-861,  1989.  . 

"Mismatch  of  patellofemoral  joint  cartilage  material  properties:  significance  m  the  etiology  of  cartilage  lesions, 
Froimson  M.I.,  Kelly  M.A.,  Gardner  T.R.,  Ratcliffe  A„  Mow  V.C.,  J.  Bone  and  Joint  Surg.,  m  review. 

"sTereopCogrammetric  determination  of  joint  anatomy  and  contact  areas  "  Soslowsky  LJ.^teshian  G.A  Mow 
V.C.,  in  Biomechanics  of  Diarthrodial  Joints,  II,  ed.  by  V.C.  Mow,  A.  Ratcliffe,  S.L.-Y.  Woo,  New  York,  Ne 

York.  Springer- Verlag,  pp.  243-268, 1990.  „ 

"Quantitation  of  articular  surface  topography  and  cartilage  thickness  in  knee  joints  using  stereophotogrammetry, 

Ateshian  G.A.,  Soslowsky  L.J.,  Mow  V.C.,  J.  Biomechanics,  24:761-776, 1991-  n  , 

"Articular  geometry  of  the  glenohumeral  joint,"  Soslowsky  L.J.,  Flatow  E.L.,  Bigliam  .  .,  ow  .  .,  in.  p. 

"Quantification  of  in  situ  contact  areas  at  the  glenohumeral  joint:  A  biomechanical  study,”  Soslowsky  L.J.,  Flatow 
E  L  Bigliani  L.U.,  Pawluk  R.J.,  Ateshian  G.A.,  Mow  V.C.,  J.  Orthop.  Res.,  10:524-534,  1992. 

"A  new  stereophotogrammetry  method  for  determining  in  situ  contact  areas  in  diarthrodial  joints:  A  comparison 
study,"  Ateshian  G.A.,  Kwak  S.D.,  Soslowsky  L.J.,  Mow  V.C.,  J.  Biomechanics  21.1 1 1-124,  1994. 

"Quantitative  anatomy  of  the  knee,"  in  Knee  Surgery,  Ateshian  G.A.,  Colman  W.W  Mow  VX.,  ed.  by  F.H.  Fu,  C.D. 

Hamer  K  G  Vince  M  D.  Miller,  Williams  and  Wilkens,  Baltimore,  MD,  pp.  101-130,  iyy4. 

"Anttiomy  and  nlchanics  of  .he  patellofemoral  join.,''  Grelsamer  R.P.,  Colmau  W.W..  Mow  V.C.,  Arthroscopy 

Review ,  2:178-188,  1994. 

Vs^DOF^nee’ testing  device  to  determine  patellar  tracking  and  patellofemoral  joint  contact  area  via 
^stereophotogrammetry,”  Gardner  T.R.,  Atesh.an  G.A.,  Grelsamer  R.P..  Mow  V.C.,  ASME  Advances  ,n 

Bioengineering  BED,  28:279-280,  1994.  .  ,  ,  .  „ v  r 

"A  new  6  degrees  of  freedom  knee  testing  apparatus  for  studies  on  experimental  surgical  techniques,  Mow  V.C., 
Henry  J.H.,  Grelsamer  R.P.,  Ahmad  C.S.,  Kwak  S.D.,  Ateshian  G.A.,  Gardner  T.R.,  presented  at  Sports  Medicine 

2000,  Stockholm,  Sweden,  June  6-8, 1995. 

"The  influence  of  iliotibial  band  tension  on  patellar  tracking  and  patellofemoral  contact,  Ahmad  C.S.,  Kwak  S.  ., 
Grelsamer  R.P.,  Henry  J.H.,  Gardner  T.R.,  Ateshian  G.A.,  Mow  V.C.,  Amer.  Acad,  of  Orthop.  Surgeons 

Annual  Meeting,  submitted. 

•^afymjS  solution  for  two  contacting  b, phasic  cartilage  layers,”  Ateshian  G.A.,  Lai  W.M.,  Zhu  W.B.,  Mow  V.C, 

J.  Biomechanics,  27:1347-1360,  1994.  .  .  .  i  .  ,  w 

"Computer  simulate  of  glenohumeral  and  patellofemoral  subluxation.Eshmatmgpa^^^ 

Flatow  E.L.,  Ateshian  G.A.,  Soslowsky  L.J.,  Pawluk  RJ.,  Grelsamer  R.P.,  Mow  V.C.,  Bigliam  L.U.,  Clin.  Orthop. 

"An  anaSmica^based  3-D  coordinate  system  for  the  knee  joint,"  Kwak  S.D  Blankevoort  L-,  Ahmad  CS  Gardner 
T.R.,  Grelsamer  R.P.,  Henry  J.H.,  Ateshian  G.A.,  Mow  V.C.,  ASME  1995  Winter  Annual  Meeting,  accepted^ 
"Development  of  multibodv  model  for  diarthrodial  joints  using  accurate  3-D  cartilage  and  bone  surfaces  Kwak  S  ., 
Ateshian  G.A.,  Blankevoort  L„  Ahmad  C.S.,  Gardner  T.R.,  Grelsamer  R.P.,  Mow  V.C.,  1995  Annual  Fall  Meeting 

of  the  Biomed.  Eng.  Society,  submitted. 
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Generator  of  Body  Data 
(GEBOD)  Program 

Huaining  Cheng 
Systems  Research  Laboratories,  Inc. 

Louise  A.  Obergefell,  Ph.D. 
Armstrong  Laboratory 

Annette  L.  Rizer 
Systems  Research  Laboratories,  Inc. 

^ a 


OBJECTIVE 


•  Automated  generation  of  data  sets  for  human  and 
dummy  rigid  body  dynamics  modeling 
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HISTORY 


•  1973  Generator  of  Occupant  Data  (GOOD) 

•  1983  Generator  of  Body  Data  (GEBOD) 

»  Use  of  anthropometric  survey  data 
»  15  segments  and  14  joints  for  male,  female,  and  child 
subjects 


•  1989  conversion  program  for  use  with  MADYMO 


HISTORY 


•  1991  GEBODHI 

»  Use  of  stereophotometric  data  for  inertial  property  and 
joint  location  calculations 
»  Separate  forearm  and  hand  segments 
»  Principal  moment  of  inertia  axes  rotated  from  local  axes 
»  Manikin  data  sets 
»  Joint  resistive  torque  properties 
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HISTORY 


•  1994  GEBODIV 

»  Correction  to  joint  center  locations 
»  New  documentation 

•  Current  version  GEBODIV.l 

»  Minor  improvements 


GEBOD  PROGRAM 


•  Interactive  program 

•  Two  methods  for  selection  subject 

»  Specify  weight  and/or  standing  height 
»  User -supplied  dimensions 

•  English  or  metric  units  for  input  and  output 

•  Separated  or  combined  hand  and  forearm  segments 
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DATA  SETS  PROVIDED 

•  Child  (2-19  years) 

•  Adult  human  male  (unit  or  percentile) 

»  Weight  range:  118  -  264  lbs 
»  Height  range:  62.17  -  77.64  in 

•  Adult  human  female  (unit  or  percentile) 

»  Weight  range:  85  -  200  lbs 
»  Height  range:  56.93  -72.05  in 


•  Hybrid  H  and  Hybrid  HI  manikins  (50th  %  tile) 


DATA  IN  DATA  SETS 


•  Segment  mass 

•  Segment  center  of  mass 

•  Segment  principal  moments  of  inertia 

•  Orientation  of  principal  axes 

•  Contact  ellipsoid  dimensions  &  location 

•  Joint  locations 

•  Joint  ranges  of  motion  &  resistive  torque 
properties 
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USER-SUPPLIED  DIMENSIONS 


No 

Dimension 

No 

Dimension 

0 

Weight 

16 

Hip  Breadth,  Standing 

1 

Standing  Height 

17 

Shoulder  to  Elbow  Length 

2 

Shoulder  Height 

18 

Forearm-Hand  Length 

3 

Armpit  Height 

19 

Biceps  Circumference 

4 

Waist  Height 

20 

Elbow  Circumference 

5 

Seated  Height 

21 

Forearm  Circumference 

6 

Head  Length 

22 

Waist  Circumference 

7 

Head  Breadth 

23 

Knee  Height,  Seated 

8 

Head  to  Chin  Height 

24 

Thigh  Circumference 

9 

Neck  Circumference 

25 

Upper  Leg  Circumference 

10 

Shoulder  Breadth 

26 

Knee  Circumference 

11 

Chest  Depth 

27 

Calf  Circumference 

12 

Chest  Breadth 

28 

Ankle  Circumference 

13 

Waist  Depth 

29 

Ankle  Height,  Outside 

14 

Waist  Breadth 

30 

Foot  Breadth 

IS 

Buttock  Depth 

31 

Foot  Length 

ADULT  OPTIONS 
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CHILD  OPTION 


10-6 


Proceedings  of  the  1995  ATB  Model  Users’  Colloquium 


STEREOPHOTOMETRIC  DATA 


•  31  male  and  46  female  subjects 

»  McConville,  et  al.,  1980;  Young,  et  al.,  1983 

•  Thousands  of  data  points 

•  Over  70  anthropometric  landmarks 

•  Segmented  by  cut  planes  defined  by  landmarks 

•  Segment  volumes  determined 

•  GEBOD  version  of  first  use:  GEBOD  in 


ANTHROPOMETRIC 
SURVEY  DATA 


•  2,420  male  and  1,905  female  subjects 

»  Grunhofer  and  Kroh,  1975;  Clauser,  et  al.,  1972 


•  3,782  children  aged  2  to  20  years 

»  Synder,  et  al.,  1977 

•  Over  150  body  dimensions  measured 

•  GEBOD  version  of  first  use:  GEBOD 
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JOINT  CENTER  LOCATION 
REGRESSION  EQUATIONS 

•  Expressions  for  joint  centers  based  on 
anthropometric  landmarks 

•  Two  separate  groups  for  male  and  female  subjects 

•  Regression  equations  are  functions  of  weight 
and/or  standing  height 


•  Used  in  GEBOD  III  to  determine  specific  joint 
center  locations  in  the  local  reference  system 


MASS  PROPERTY 
REGRESSION  EQUATIONS 


•  Based  on  stereophotometric  segment  volumes 

•  Two  separate  groups  for  male  and  female  subjects 

•  Functions  of  weight  and/or  standing  height 

•  Used  in  GEBOD  HI  to  predict  each  segment’s  mass 
and  principal  moments  of  inertia 
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BODY  DIMENSION  SET 
REGRESSION  EQUATIONS 


•  Based  on  anthropometric  survey 

•  Two  separate  groups  for  male  and  female 

•  Functions  of  weight  and/or  standing  height 

»  Extra  independent  variable  of  age  for  child  regression  equations 


BODY  DIMENSION  SET 
REGRESSION  EQUATIONS 


•  32  body  dimensions  used 

»  Three  additional  dimensions  used  when  separate  hands  modeled 
#32:  Hand  breadth;  #33:  Hand  length;  #34  Hand  depth 

•  Used  in  GEBODIII  to  determine  contact  ellipsoid 
dimensions  and  locations 
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EXAMPLE  BODY  CONFIGURATION 


Head  Pivot 
Neck  Pivot 

Left/Right  Shoulder 
Waist 

Left/Right  Elbow 
Pelvis 
Left/Right  Hip 
Left/Right  Wrist 


Left/Right  Knee 


Left/Right  Ankle 


Head 

Neck 

Upper  Torso 
LeR/Right  Upper  Arm 
Middle  Torso 

Left/Right  Lower  Arm 
Lower  Torso 
Left/Right  Hand 
Left/Right  Upper  Leg 


Left/Right  Lower  Leg 
Left/Right  Foot 


JOINT  RESISTIVE 
CHARACTERISTICS 


•  Resistive  torque  properties  for  spring  torque, 
viscous  torque,  and  coulomb  friction 

•  Joint  lock/unlock  conditions 

•  Joint  stop  angles 

•  Obtained  from  a  number  of  sources 
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OUTPUT  FILES 


•  Extension  .TAB  file 

»  A  list  of  35  body  dimensions  for  the  human  data  sets 
»  Labeled  tables  of  the  calculated  data  set 

•  Extension  .AIN  file 

»  ATB  model  input  cards  B.l  to  B.6  in  input  ready  format 


SUMMARY 


•  Child,  adult  male  and  female,  and  Hybrid  dummy 
data  sets 

•  Interactive  and  flexible  in  amount  of  user  input 


•  Provides  the  body  data  necessary  for  human  body 
dynamic  simulations 
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C : \ATB>gebodivl 


PROGRAM  GEBOD 

GEBOD  GENERATES  BODY  DESCRIPTION  DATA 
SUITABLE  FOR  INPUT  TO  THE  ATB  MODEL 

PLEASE  ENTER  A  DESCRIPTION  OF  THE  SUBJECT  ( <60  CHARS.) 
male  72  inches  tall  and  195  lbs 


ENTER  FILENAME  FOR  ALL  OUTPUT  FILES, 

EXTENSIONS  WILL  BE  ASSIGNED: 
ucp\m72-195 

1)  CHILD  (2-19  YEARS) 

2)  ADULT  HUMAN  FEMALE 

3)  ADULT  HUMAN  MALE 

4)  USER-SUPPLIED  BODY  DIMENSIONS 

5)  SITTING  HYBRID  III  DUMMY  (50%) 

6)  STANDING  HYBRID  III  DUMMY  (50%) 

7)  HYBRID  II  DUMMY  (50%) 

ENTER  NUMBER  CORRESPONDING  TO  DESIRED  SUBJECT  TYPE 
3 

1)  WEIGHT 

2)  STANDING  HEIGHT 

3)  ALL  OF.  THE  ABOVE 

ENTER  NUMBER  CORRESPONDING  TO 
PREDICTING  DIMENSION ( S )  TO  BE  SUPPLIED 
3 


SELECT  UNITS  FOR  WEIGHT 

1)  LB. 

2)  N. 

3)  %-TILE 

1 

ENTER  VALUE  FOR  WEIGHT 

IN  THE  RANGE  118.0  ,  264.0 

195 
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SELECT  UNITS  FOR  STANDING  HEIGHT 

1)  IN. 

2)  M. 

3)  %-TILE 

1 

ENTER  VALUE  FOR  STANDING  HEIGHT 
IN  THE  RANGE  62.17  ,  77.64 

72 


1)  FOREARM  AND  HAND  COMBINED 

2 )  FOREARM  AND  HAND  SEPARATED 

enter  the  number  corresponding  to  the  desired  lower  arm  segmentation 
1 


SELECT  UNITS  FOR  OUTPUT 

1)  ENGLSH 

2)  METRIC 

1 


IS  ATB  MODEL  FORMATTED  OUTPUT  DESIRED  (Y/N)  ? 

y 


IS  IT  DESIRED  TO  PRODUCE  ANOTHER  BODY  DESCRIPTION  DATA  SET  (Y/N)  ? 
n 

FROM  ARESULTS 
C : \ATB> 
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C: \ATB\UCP>ndeckl 

THE  DEFAULT  RESPONSE  TO  QUESTIONS  IN  THIS  PROGRAM  IS 
<YES> .  A  CR  (CARRIAGE  RETURN ) ,  Y /  OR  YES  RESPONSE 
WILL  BE  INTERPRETED  AS  AN  AFFIRMATIVE  RESPONSE. 

ANY  OTHER  RESPONSE  WILL  BE  INTERPRETED  AS  A 
NEGATIVE  RESPONSE. 

ENTER  COMMENT  80  MAX.  PER  LINE  160  CHRS  IN  TWO  LINES 

Example  setup  file  using  a  standard  vehicle  interior  and  static  pulse 
Occupant  from  GEBOD,  male  72  inches  and  195  lbs 
ENTER  TOTAL  SIMULATION  TIME  (SEC) 

.2 


1.  FIRST  COMMENT  - 

Example  setup  file  using  a  standard  vehicle  interior  and  static  pulse 

2.  SECOND  COMMENT  - 

Occupant  from  GEBOD,  male  72  inches  and  195  lbs 

3.  TOTAL  SIMULATION  TIME  =  .200 

4.  NPRT  ARRAY  - 

12345678  9101112131415161718192021222324252627282930313233343536 
1  010  20000  00  00000000000000000000000001 
THESE  ARE  DEFAULT  NPRT  VALUES 

DO  YOU  WISH  TO  MAKE  ANY  CHANGES  ?  <YES> 
n 


DIRECTORY  OF  OCCUPANT  DATA  LIBRARY 


FILE  NO.  FILE  NAME 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


SIXYRCHD 

95THMALE 

PART572 

ESTP572 

TEST  50-50-M 

PED50TH 

75-75-M 

25-25-M 

5-5-F 

95THMLOCC 

50-50-M 

15-95-M 

37-37-F 


ENTER  THE  DESIRED  FILE  NO.  FOR  INPUT  DATA 
ENTER  0  TO  USE  A  GEBOD  FILE 
0 

ENTER  THE  FILENAME  OF  THE  GEBOD  FILE  TO  USE 
m72-195 . ain 
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DIRECTORY  OF  VEHICLE  DECELERATION  DATA  LIBRARY 


FILE  NO.  FILE  NAME 

OG  STATIC 
RX-7  NCAP 
PULSE  #1 
OMNI  SLED 
VOLVO  BARRIER 
SIDE  IMP  (DRIVER 
FIAT  DATA 
SUBARU  BRAT 
G30F 


ENTER  THE  DESIRED  FILE  NO.  FOR  INPUT  DATA 

VEHICLE  DECELERATION  CHOSEN  IS  : 

OG  STATIC 

IS  THIS  O.K.  ?  <YES> 

Y  _ 

DIRECTORY  OF  VEHICLE  INTERIOR  DATA  LIBRARY 

FILE  NO.  FILE  NAME 

MAZDA  RX-7 
'88  CHEV.  BERETT 
STD  VEH  #  1 
VOLVO 
OMNI -MOD 
OMNI-STD 
OMNI -STD ( UNR ) 

SIDE  IMPACT  BASE 
PEDVEH1 
STD  VEH  #  2 
CONTACT  CHK 
G24 

FORKLIFT 


ENTER  THE  DESIRED  FILE  NO.  FOR  INPUT  DATA 
10 

VEHICLE  INTERIOR  CHOSEN  IS  : 

STD  VEH  #  2 

IS  THIS  O.K.  ?  < YES> 


Stop  -  Program  terminated. 


C : \ATB\UCP> 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 


1 

2 

3 

4 

5 

6 

7 

8 
9 
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C: \ATB\UCP>ndeck2 

DO  YOU  WANT  ALL  SEGMENTS  TO  BE  RIGID  (YES  OR  NO )  ?  : 
n 

WILL  ALL  SEGMENTS  HAVE  THE  SAME  PROPERTIES  ?: 

Y  _ _ 

*****  ALL  SEGMENTS  WILL  HAVE  THIS  VALUE  ***** 
FUNCTION  SELECTION  -  SEGMENT  NO.  1  ,  LT 
FORCE  -  DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER  (999  TO  SEE  LIST): 

999 


FUNC.  NO.  FUNC.  NAME 

1  SEAT  BOTTOM  FDF 

2  SEAT  BACK  FDF 

3  STIFF  SURFACES  FDF 

4  DASH  -  FEMUR  FDF 

5  FRONT  DASH-FEMUR  FDF 

6  DASH  -  TORSO  FDF 

7  HEAD  -  WINDSHLD  FDF 

8  COLUMN  -  BODY  FDF 

9  STEERING  COL  E/A  FDF 

10  SEGMENT-SEGMENT  FDF 

11  HARNESS  (CABLE)  FDF 

12  HARNESS  (LAP)  FDF 

13  HARNESS  ( SHLDR )  FDF 

14  HARNESS  FDF 
FORCE  -  DEFLECTION  FUNCTION 

ENTER  FUNCTION  NUMBER  (999  TO  SEE  LIST): 


*****  ALL  SEGMENTS  WILL  HAVE  THIS  VALUE  ***** 
FUNCTION  SELECTION  -  SEGMENT  NO.  1  ,  LT 
ENERGY  ABSORPTION  (R)  FUNCTION 
ENTER  999  TO  SEE  LIST  OR  FUNCTION  NUMBER  : 

999 


FUNC.  NO.  FUNC.  NAME 

1  SEAT  BOTTOM  R=0.625 

2  SEAT  BACK  R  =  0.53 

3  STIFF  SURF.  R  =  1.0 

4  DASH  -  FEMUR  R=0.124 

5  DASH  -  TORSO  R=0.133 

6  STEERING  COLMN  R=0 . 0 

7  SEG.  -  SEG.  R- FACTOR 
ENERGY  ABSORPTION  (R)  FUNCTION 

ENTER  999  TO  SEE  LIST  OR  FUNCTION  NUMBER  : 
7 


*****  ALL  SEGMENTS  WILL  HAVE  THIS  VALUE  ***** 
FUNCTION  SELECTION  -  SEGMENT  NO.  1  ,  LT 
PERMANENT  SET  (G)  FUNCTION 

ENTER  999  TO  SEE  LIST  OR  FUNCTION  NUMBER  : 

999 
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FUNC .  NO . 


FUNC .  NAME 


1  SEAT  G  =  0.08 

2  STIFF  SURF.  G  =  0.88 

3  DASH  -  FEMUR  G=0.809 

4  DASH  -  TORSO  G=0.597 

5  SEG-SEG  CONTCT  G=0 . 0 
PFRMANENT  SET  (G)  FUNCTION 

ENTER  999  TO  SEE  LIST  OR  FUNCTION  NUMBER 

5 


FUNC.  NO.  FUNC.  NAME 

1  SEAT  BOTTOM  FDF 

2  SEAT  BACK  FDF 

3  STIFF  SURFACES  FDF 

4  DASH  -  FEMUR  FDF 

5  FRONT  DASH-FEMUR  FDF 

6  DASH  -  TORSO  FDF 

7  HEAD  -  WINDSHLD  FDF 

8  COLUMN  -  BODY  FDF 

9  STEERING  COL  E/A  FDF 

10  SEGMENT-SEGMENT  FDF 

11  HARNESS  (CABLE)  FDF 

12  HARNESS  (LAP)  FDF 

13  HARNESS  ( SHLDR )  FDF 

14  HARNESS  FDF 


FUNCTION  SELECTION  BY  SEGMENT 
SEGMENT  FUNCTION  NO. 


DO  YOU  WANT  TO  MAKE  ANY  CHANGES  (YES  OR^NO, 
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FUNC .  NO . 


FUNC.  NAME 


SEAT  BOTTOM  R=0.625 
SEAT  BACK  R  =  0.53 
STIFF  SURF.  R  =  1.0 
DASH  -  FEMUR  R=0.124 
DASH  -  TORSO  R=0.133 
STEERING  COLMN  R=0.0 
SEG.  -  SEG.  R-FACTOR 


FUNCTION  SELECTION  BY  SEGMENT 
SEGMENT  FUNCTION  NO. 

1  LT  7 


11  LF 

12  RUA 


DO  YOU  WANT  TO  MAKE  ANY  CHANGES  (YES  OR  NO)  n 


FUNC.  NO.  FUNC.  NAME 


SEAT  G  =  0.08 
STIFF  SURF.  G  =  0.88 
DASH  -  FEMUR  G=0.809 
DASH  -  TORSO  G= 0.597 
SEG-SEG  CONTCT  G=0.0 


FUNCTION  SELECTION  BY  SEGMENT 
SEGMENT  FUNCTION  NO. 


DO  YOU  WANT  TO  MAKE  ANY  CHANGES  (YES  OR  NO )  n 
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FUNCTION  SELECTION 

PLANE  NO.  1  SEAT  CUSHION 


FORCE-DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  0  FOR  NO  CONTACTS, 

1 


999  TO  SEE  LIST)  : 


FUNCTION  SELECTION 
PLANE  NO.  1  SEAT  CUSHION 
INERTIAL  SPIKE  FUNCTION 
ENTER  FUNCTION  NUMBER 
ENTER  0  FOR  NO  INERTIAL  SPIKE, 
0 


999  TO  SEE  LIST) 


FUNCTION  SELECTION 

PLANE  NO.  1  SEAT  CUSHION 

ENERGY  DISSIPATION  (R)  FUNCTION 


ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 


FUNCTION  SELECTION 
PLANE  NO.  1  SEAT  CUSHION 

PERMANENT  SET  (G)  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

1 

FUNCTION  SELECTION 
PLANE  NO.  1  SEAT  CUSHION 

FRICTION  COEFFICIENT  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

2 


FUNCTION  SELECTION 
PLANE  NO.  2  SEAT  BACK 


FORCE-DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  0  FOR  NO  CONTACTS, 
2 


999  TO  SEE  LIST)  : 


FUNCTION  SELECTION 
PLANE  NO.  2  SEAT  BACK 
INERTIAL  SPIKE  FUNCTION 
ENTER  FUNCTION  NUMBER 
ENTER  0  FOR  NO  INERTIAL  SPIKE, 
0 


999  TO  SEE  LIST) 


FUNCTION  SELECTION 
PLANE  NO.  2  SEAT  BACK 
ENERGY  DISSIPATION  (R)  FUNCTION 


ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 


FUNCTION  SELECTION 
PLANE  NO.  2  SEAT  BACK 
PERMANENT  SET  (G)  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

lf-8 
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FUNCTION  SELECTION 
PLANE  NO.  2  SEAT  BACK 
FRICTION  COEFFICIENT  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

2 


FUNCTION  SELECTION 

PLANE  NO.  3  FLOOR  BOARD 


FORCE-DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 

(ENTER  0  FOR  NO  CONTACTS,  999  TO  SEE  LIST)  : 
3 


FUNCTION  SELECTION 
PLANE  NO.  3  FLOOR  BOARD 
INERTIAL  SPIKE  FUNCTION 
ENTER  FUNCTION  NUMBER 

ENTER  0  FOR  NO  INERTIAL  SPIKE,  999  TO  SEE  LIST)  : 
0 


FUNCTION  SELECTION 
PLANE  NO.  3  FLOOR  BOARD 
ENERGY  DISSIPATION  (R)  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

3 


FUNCTION  SELECTION 
PLANE  NO.  3  FLOOR  BOARD 
PERMANENT  SET  (G)  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

2 


FUNCTION  SELECTION 
PLANE  NO.  3  FLOOR  BOARD 
I  FRICTION  COEFFICIENT  FUNCTION 

ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

3 


FUNCTION  SELECTION 
PLANE  NO.  4  TOE  BOARD 


FORCE-DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 

(ENTER  0  FOR  NO  CONTACTS,  999  TO  SEE  LIST)  : 
3 


FUNCTION  SELECTION 
PLANE  NO.  4  TOE  BOARD 
INERTIAL  SPIKE  FUNCTION 
ENTER  FUNCTION  NUMBER 

ENTER  0  FOR  NO  INERTIAL  SPIKE,  999  TO  SEE  LIST)  : 
0 
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FUNCTION  SELECTION 
PLANE  NO.  4  TOE  BOARD 
ENERGY  DISSIPATION  (R)  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

3 


FUNCTION  SELECTION 
PLANE  NO.  4  TOE  BOARD 
PERMANENT  SET  (G)  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

2 


FUNCTION  SELECTION 
PLANE  NO.  4  TOE  BOARD 
FRICTION  COEFFICIENT  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  : 

3 


FUNCTION  SELECTION 
PLANE  NO.  5  WINDSHIELD 


FORCE-DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 

(ENTER  0  FOR  NO  CONTACTS,  999  TO  SEE  LIST)  : 
0 

PLANE  SELECTED  AS  HAVING  NO  CONTACTS 


FUNCTION  SELECTION 
PLANE  NO.  6  ROOF 


FORCE-DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  0  FOR  NO  CONTACTS, 


999  TO  SEE  LIST) 


0 

PLANE  SELECTED  AS 


HAVING  NO  CONTACTS 


FUNCTION  SELECTION 
PLANE  NO.  7  DASH 


FORCE-DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  0  FOR  NO  CONTACTS,  999 


TO  SEE  LIST) 


0 

PLANE  SELECTED  AS  HAVING  NO 


CONTACTS 


FUNCTION  SELECTION 

PLANE  NO.  8  KNEE  BOLSTER 


FORCE-DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  0  FOR  NO  CONTACTS, 


999  TO  SEE  LIST)  : 


0 

PLANE  SELECTED  AS  HAVING  NO 


CONTACTS 
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FUNC.  NO.  FUNC.  NAME 

1  SEAT  BOTTOM  FDF 

2  SEAT  BACK  FDF 

3  STIFF  SURFACES  FDF 

4  DASH  -  FEMUR  FDF 

5  FRONT  DASH-FEMUR  FDF 

6  DASH  -  TORSO  FDF 

7  HEAD  -  WINDSHLD  FDF 

8  COLUMN  -  BODY  FDF 

9  STEERING  COL  E/A  FDF 

10  SEGMENT-SEGMENT  FDF 

11  HARNESS  (CABLE)  FDF 

12  HARNESS  (LAP)  FDF 

13  HARNESS  ( SHLDR )  FDF 

14  HARNESS  FDF 

FUNCTION  SELECTION  BY  PLANE 
PLANE  FUNCTION  NO. 

1  SEAT  CUSHION  1 

2  SEAT  BACK  2 

3  FLOOR  BOARD  3 

4  TOE  BOARD  3 

5  WINDSHIELD 

6  ROOF 

7  DASH 

8  KNEE  BOLSTER 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES  (YES  OR  NO )  n 


FUNC.  NO.  FUNC.  NAME 

1  HEAD  -  WINDSHLD  SPIK 

FUNCTION  SELECTION  BY  PLANE 
PLANE  FUNCTION  NO. 

1  SEAT  CUSHION 

2  SEAT  BACK 

3  FLOOR  BOARD 

4  TOE  BOARD 

5  WINDSHIELD 

6  ROOF 

7  DASH 

8  KNEE  BOLSTER 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES  (YES  OR  NO )  n 


FUNC.  NO.  FUNC.  NAME 

1  SEAT  BOTTOM  R=0.625 

2  SEAT  BACK  R  =  0.53 

3  STIFF  SURF.  R  =  1.0 

4  DASH  -  FEMUR  R=0.124 

5  DASH  -  TORSO  R=0.133 

6  STEERING  COLMN  R=0 . 0 

7  SEG.  -  SEG.  R-FACTOR 

FUNCTION  SELECTION  BY  PLANE 
PLANE  FUNCTION  NO. 

1  SEAT  CUSHION 

2  SEAT  BACK 

3  FLOOR  BOARD 


4  TOE  BOARD  3 

5  WINDSHIELD  0 

6  ROOF  0 

7  DASH  0 


8  KNEE  BOLSTER  n  n  0 

Pro^g^&U iuANY  CHANGES  (YES  OR  NO)  n 


OOOOI  OOOOOOOOI  H  (N  m 


FUNC.  NO. 


FUNC .  NAME 


1  SEAT  G  =  0.08 

2  STIFF  SURF.  G  =  0.88 

3  DASH  -  FEMUR  G=0.809 

4  DASH  -  TORSO  G=0.597 

5  SEG-SEG  CONTCT  G=0.0 


DO  YOU  WANT  TO  MAKE  ANY  CHANGES 


1 

2 

3 

4 

5 

6 

7 

8 

(YES  OR  NO) 


FUNCTION  SELECTION  BY  PLANE 
PLANE  FUNCTION  NO. 

SEAT  CUSHION 
SEAT  BACK 
FLOOR  BOARD 
TOE  BOARD 
WINDSHIELD 
ROOF 
DASH 

KNEE  BOLSTER 
n 


1 

1 

2 

2 

0 

0 

0 

0 


FUNC. 

NO .  FUNC .  NAME 

1 

FRICTION  MU=1 . 0 

2 

FRICTION  MU  =  0.75 

3 

FRICTION  MU=0 . 5 

4 

HARNESS  MU  =  0.2 

5 

HARNESS  MU  =  0.1 

6 

CONSTANT,  F  =  0 . 0 

FUNCTION  SELECTION  BY  PLANE 

PLANE  FUNCTION  NO. 

1 

SEAT  CUSHION 

2 

SEAT  BACK 

3 

FLOOR  BOARD 

4 

TOE  BOARD 

5 

WINDSHIELD 

6 

ROOF 

7 

DASH 

8 

KNEE  BOLSTER 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES 

(YES  OR  NO) 

n 
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(NCNCOCQOOOO 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE  1 
SEAT  CUSHION  : 

3 

ENTER  SEGMENT  NOS.  FOR  CONTACT  WITH  THIS  PLANE 

CONTACT  NO .  1 

INPUT  SEGMENT  NO .  :  1 

CONTACT  NO .  2 

INPUT  SEGMENT  NO .  :  6 

CONTACT  NO .  3 

INPUT  SEGMENT  NO.  :  9 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE  2 
SEAT  BACK  : 

3 

ENTER  SEGMENT  NOS.  FOR  CONTACT  WITH  THIS  PLANE 

CONTACT  NO .  1 

INPUT  SEGMENT  NO.  :  1 

CONTACT  NO .  2 

INPUT  SEGMENT  NO .  :  2 

CONTACT  NO .  3 

INPUT  SEGMENT  NO .  :  3 
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SEGMENT  NO.  TITLE 

1  LT 

2  CT 

3  UT 

4  N 

5  H 

6  RUL 

7  RLL 

8  RF 

9  LUL 

10  LLL 

11  LF 

12  RUA 

13  RLA 

14  LUA 

15  LLA 


PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE  3 
FLOOR  BOARD  : 


2 

ENTER  SEGMENT  NOS 
CONTACT  NO .  1 

INPUT  SEGMENT  NO. 
CONTACT  NO .  2 

INPUT  SEGMENT  NO. 


FOR  CONTACT  WITH  THIS  PLANE 
:  8 
:  11 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE  4 
TOE  BOARD  : 

2 

ENTER  SEGMENT  NOS.  FOR  CONTACT  WITH  THIS  PLANE 

CONTACT  NO.  1 

INPUT  SEGMENT  NO.  :  8 

CONTACT  NO.  2 

INPUT  SEGMENT  NO.  :  11 
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PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE 
WINDSHIELD  = 

0 

NO  CONTACTS  ALLOWED  FOR  THIS  PLANE 


5 


PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE  6 
ROOF  : 

0 

NO  CONTACTS  ALLOWED  FOR  THIS  PLANE 


PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE 
DASH  : 
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NO  CONTACTS  ALLOWED  FOR  THIS  PLANE 


SEGMENT  NO.  TITLE 

1  LT 

2  CT 

3  UT 

4  N 

5  H 

6  RUL 

7  RLL 

8  RF 

9  LUL 

10  LLL 

11  LF 

12  RUA 

13  RLA 

14  LUA 

1 5  LLA 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  NO.  TITLE 


1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

11-15 

ment  no. 

title 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE  8 
KNEE  BOLSTER  : 

0 

NO  CONTACTS  ALLOWED  FOR  THIS  PLANE 


PLANE 

NO.  TITLE 


SEGMENT  CONTACTS 
SEGMENTS 


1  SEAT  CUSHION  3:  1-  LT  ,  6-  RUL ,  , 

2  SEAT  BACK  3:  1~  LT  '  I  '  3  ' 

3  FLOOR  BOARD  2:  8-  RF  ,11  LF  , 

4  TOE  BOARD  2:  8-  RF  ,  11-  LF  , 

5  WINDSHIELD  0: 

6  ROOF  0 : 

7  DASH  0 : 

8  KNEE  BOLSTER  0:  „„ .  „ 

DO  YOU  WISH  TO  MODIFY  ANY  CONTACTS  (YES  OR  NO )  . 

DO  YOU  WISH  TO  SPECIFY  ANY  SEGMENT-SEGMENT  CONTACTS  (YES  OR  NO 
Y  _ 


9-  LUL, 
3-  UT  , 


8-  RF 
8-  RF 


11-  LF 
11-  LF 


MENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION  , 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  1 

LT  : 

0 
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SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  2 
CT  : 

0 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  3 
UT  : 

0 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  4 
N  : 

0 
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InterNnoCOofAsegmentsTcontacting  SEGMENT  6 
RUL  : 

ENTER  SEGMENT  NOS.  FOR  CONTACT  WITH  THIS  SEGMENT 


CONTACT  NO .  1 

INPUT  SEGMENT  NO.  :  9 


CONTACT  NO .  2 

INPUT  SEGMENT  NO.  :  13 
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SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  7 
RLL  : 

2 

ENTER  SEGMENT  NOS.  FOR  CONTACT  WITH  THIS  SEGMENT 

CONTACT  NO .  1 

INPUT  SEGMENT  NO.  :  10 

CONTACT  NO .  2 

INPUT  SEGMENT  NO.  :  11 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

.  RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  8 
RF  : 

2 

ENTER  SEGMENT  NOS.  FOR  CONTACT  WITH  THIS  SEGMENT 

CONTACT  NO.  1 

INPUT  SEGMENT  NO.  :  10 

CONTACT  NO.  2 

INPUT  SEGMENT  NO.  :  11 
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SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  9 
LUL  : 

ENTER  SEGMENT  NOS.  FOR  CONTACT  WITH  THIS  SEGMENT 

CONTACT  NO .  1 

INPUT  SEGMENT  NO.  :  15 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  10 
LLL  : 

0 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 
ENTER  NO.  OF  SEGMENTS  CONTACTING 
LF  : 

0 
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SEGMENT  NO.  TITLE 

1  LT 

2  CT 

3  UT 

4  N 

5  H 

6  RUL 

7  RLL 

8  RF 

9  LUL 

1 0  LLL 

11  LF 

12  RUA 

13  RLA 

14  LUA 

15  LLA 


SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  12 
RUA  : 

0 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  13 
RLA  : 

0 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  14 
LUA  : 

0 
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SEGMENT  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  SEGMENT  15 
LLA  : 

0 


SEGMENT  NO. 

TITLE 

1 

LT 

2 

CT 

3 

UT 

4 

N 

5 

H 

6 

RUL 

7 

RLL 

8 

RF 

9 

LUL 

10 

LLL 

11 

LF 

12 

RUA 

13 

RLA 

14 

LUA 

15 

LLA 

SEGMENT 
NO.  TITLE 


SEGMENT  CONTACTS 
NO.  SEGMENTS 


8  RF 

9  LUL 

10  LLL 

11  LF 

12  RUA 

13  RLA 

14  LUA 

15  LLA 


9-  LUL,  13-  RLA, 
10-  LLL,  11-  LF  , 
10-  LLL,  11-  LF  , 
15-  LLA, 


DO  YOU  WISH  TO  MODIFY  ANY  CONTACTS  (YES  OR  NO )  ? 


FUNCTION  SELECTION  -  SEGMENT  NO.  6  ,  RUL 

FRICTION  COEFFICIENT  (C)  FUNCTION 

ENTER  FUNCTION  NUMBER  (999  TO  SEE  LIST): 

3 


FUNCTION  SELECTION  -  SEGMENT  NO.  7  ,  RLL 

FRICTION  COEFFICIENT  (C)  FUNCTION 

ENTER  FUNCTION  NUMBER  (999  TO  SEE  LIST): 

3 


FUNCTION  SELECTION  -  SEGMENT  NO.  8  ,  RF 

FRICTION  COEFFICIENT  (C)  FUNCTION 

ENTER  FUNCTION  NUMBER  (999  TO  SEE  LIST): 

3 


FUNCTION  SELECTION  -  SEGMENT  NO.  9  ,  LUL 

FRICTION  COEFFICIENT  (C)  FUNCTION 

ENTER  FUNCTION  NUMBER  (999  TO  SEE  LIST): 

3 
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FUNC.  NO. 


FUNC.  NAME 


1  FRICTION  MU=1 . 0 

2  FRICTION  MU  =  0.75 

3  FRICTION  MU=0 . 5 

4  HARNESS  MU  =  0.2 

5  HARNESS  MU  =  0.1 

6  CONSTANT,  F=0.0 

FUNCTION  SELECTION  BY  SEGMENT 
SEGMENT  FUNCTION  NO. 

1  LT 

2  CT 

3  UT 

4  N 

5  H 

6  RUL 

7  RLL 

8  RF 

9  LUL 

10  LLL 

11  LF 

12  RUA 

13  RLA 

14  LUA 

15  LLA 

DO  YOU  WANT  TO  MAKE  ANY  CHANGES  (YES  OR  NO )  n 


Stop  -  Program  terminated. 


C : \ATB\UCP> 
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ooooocomcocooooooo 


C • \ATB\UCP>ndeck3 

THE  DEFAULT  RESPONSE  TO  QUESTIONS  IN  THIS  PROGRAM  IS 
<YES> .  A  <CR>,  Y,  OR  YES  RESPONSE'  WILL  BE 
INTERPRETED  AS  AN  AFFIRMATIVE  RESPONSE. 

FOR  A  NEGATIVE  RESPONSE,  YOU  MUST  ENTER  "N"  OR  "n" 


INPUT  FOR  THE  DEFAULT  OPTION  HAS  BEEN  GENERATED 
DO  YOU  ACCEPT  THE  DEFAULT  OPTION?  <YES> 

INPUT  DATA  WILL  BE  LISTED  FOR  ENTRY  CORRECTION 
IS  THIS  O.K.  ?  <YES> 

Y 


SEGMENT  LINEAR 
SEGMENT  NO. 

5 

3 

1 

8 

11 

16 


ACCELERATIONS 


000 

LOCATION (X,Y,Z) 
.000 

.000 

000 

.000 

.  000 

000 

.000 

.000 

000 

.000 

.000 

000 

.000 

.000 

000 

.000 

.  000 

DO  YOU  ACCEPT  THESE  INPUTS  ?  <YES> 

IF  NOT,  THE  INPUTS  MAY  BE  RESPECIFIED 


SEGMENT  LINEAR  VELOCITIES 

SEGMENT  NO.  LOCATION ( X , Y , Z ) 

5  .000  .000  .000 

3  .000  .000  .000 

1  .000  .000  .000 

8  .000  .000  .000 

DO  YOU  ACCEPT  THESE  INPUTS  ?  <YES> 

IF  NOT,  THE  INPUTS  MAY  BE  RESPECIFIED 

y 


SEGMENT  LINEAR  DISPLACEMENTS 

SEGMENT  NO.  LOCATION ( X , Y , Z ) 


000 

.  000 

000 

.000 

000 

.000 

000 

.000 

000 

.000 

000 

.000 

DO  YOU  ACCEPT  THESE  INPUTS  ?  <YES> 

IF  NOT,  THE  INPUTS  MAY  BE  RESPECIFIED 

y 


.000 
.000 
.000 
.  000 

.  000 

.000 


SEGMENT  ANGULAR  ACCELERATIONS 
SEGMENT  NO. 

5 

3 

1 

DO  YOU  ACCEPT  THESE  INPUTS  ?  <YES> 

IF  NOT,  THE  INPUTS  MAY  BE  RESPECIFIED 

y 
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SEGMENT  ANGULAR  VELOCITIES 
SEGMENT  NO. 

5 

3 
1 

DO  YOU  ACCEPT  THESE  INPUTS  ?  <YES> 

IF  NOT,  THE  INPUTS  MAY  BE  RESPECIFIED 

y 

SEGMENT  ANGULAR  DISPLACEMENTS 
SEGMENT  NO. 

5 

4 
3 
2 
1 

7 
9 

8 

10 

DO  YOU  ACCEPT  THESE  INPUTS  ?  <YES> 

IF  NOT,  THE  INPUTS  MAY  BE  RESPECIFIED 

y 

JOINT  PARAMETERS 
JOINT  NO. 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

DO  YOU  ACCEPT  THESE  INPUTS  ?  <YES> 

IF  NOT,  THE  INPUTS  MAY  BE  RESPECIFIED 

Y 

Stop  -  Program  terminated. 


C : \ATB\UCP> 
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C :  \ATB\UCP>NDECK4  to 

THE  DEFAULT  RESPONSE  TO  QUESTIONS  IN  THIS  PROGRAM  IS 
<YES>  .  A  <CR> ,  Y ,  OR  YES  RESPONSE'  WILL  BE 
INTERPRETED  AS  AN  AFFIRMATIVE  RESPONSE.  ANY  OTHER 
RESPONSE  WILL  BE  INTERPRETED  AS  A  NEGATIVE  RESPONSE. 

PLEASE  WAIT  -  LOADING  DATA  FILES 


CLEARANCE  BETWEEN  FOOT  AND  TOEBOARD  IS  ONE  INCH 
DO  YOU  WISH  TO  CHANGE  THIS  VALUE  ?<YES> 

N 

DO  YOU  ACCEPT  THIS  INITIAL  CONDITION  ? 


N 

FOLLOWING  OPTIONS  ARE  AVAILABLE 


ENTER  D  FOR  SEGMENT  POSITION  AND 
E  FOR  ENTERING  NEW  POSITION 
Q  FOR  QUITTING 


ORIENTATIONS 
AND  ORIENTATIONS 


D 
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SEG 

X 

-6.22 
-7.31 
-12 . 08 
-14.94 
-13.67 

1.55 
16.45 

27.55 

1.55 
16.45 

27.55 
-13.21 
-4.35 
-13.21 
-4.35 


SEG  ORIENTATION 
PITCH  ROLL 

00  25.24  .00 


Y 

Z 

YAW 

PITCH 

.00 

-4.93 

.00 

25.24 

.00 

-9.50 

.00 

25.24 

.  00 

-16.06 

.00 

25.24 

.00 

-24.31 

.00 

12.73 

.00 

-28.92 

.00 

.  00 

4.31 

-5.61 

.  00 

106.78 

5.37 

-3.29 

.00 

50.69 

4.95 

1.58 

.00 

140.69 

-4.31 

-5.61 

.00 

106.78 

-5.37 

-3.29 

.00 

50.69 

-4.95 

1.58 

.00 

140.69 

6.98 

-14.70 

.00 

25.02 

6.21 

-10.30 

.00 

90.00 

-6.98 

-14.70 

.00 

25.02 

-6.21 

-10.30 

.00 

90.00 

:ONS  ARE 

AVAILABLE 

ENTER  D  FOR  SEGMENT  POSITION  AND  ORIENTATIONS 

E  FOR  ENTERING  NEW  POSITION  AND  ORIENTATIONS 
Q  FOR  QUITTING 
E 

ENTER  SEGMENT  NO.  TO  BE  CHANGED,  TO  STOP,  ENTER  0 
13 

ENTER  VARIABLE  TYPE  :  1.  LINEAR 

2.  ANGULAR 

2 

ENTER  THE  VARIABLE  TYPE  INDEX 
^1,2,3  FOR  X , Y , Z  OR  YAW, PITCH , ROLL ) 

ENTER  THE  NEW  VALUE  OF  THE  VARIABLE 
105 

ENTER  SEGMENT  NO.  TO  BE  CHANGED,  TO  STOP,  ENTER  0 
15 

ENTER  VARIABLE  TYPE  :  1.  LINEAR 

2 .  ANGULAR 

2 

ENTER  THE  VARIABLE  TYPE  INDEX 
(1,2,3  FOR  X,Y,Z  OR  YAW, PITCH , ROLL ) 

2 

ENTER  THE  NEW  VALUE  OF  THE  VARIABLE 
105 

ENTER  SEGMENT  NO.  TO  BE  CHANGED,  TO  STOP,  ENTER  0 
0 

FOLLOWING  OPTIONS  ARE  AVAILABLE 

ENTER  D  FOR  SEGMENT  POSITION  AND  ORIENTATIONS 

E  FOR  ENTERING  NEW  POSITION  AND  ORIENTATIONS 
Q  FOR  QUITTING 
D 


SEG 

X 

-6.22 
-7.31 
-12.08 
-14.94 
-13.67 
1.55 
16.45 
27.55 
1 .55 
16.45 
27 . 55 
-13.21 
-4 . 74 
-13.21 
-4.74 


POSITION 

Y 

.00 

.00 

.00  - 
.00  - 
.00  - 
4.31 
5.37 
4.95 
-4.31 
-5.37 
-4.95 
6.98  - 

6.21  - 

-6.98  - 

-6.21  - 


Z 

-4.93 

-9.50 

-16.06 

-24.31 

-28.92 

-5.61 

-3.29 

1.58 

-5.61 

-3.29 

1.58 

-14.70 

-12.16 

-14.70 

-12.16 


SEG 

YAW 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 

.00 


ORIENTATION 


PITCH 
25.24 
25.24 
25.24 
12.73 
.  00 
106.78 
50.69 
140 . 69 
106.78 
50.69 
140.69 
25.02 
105.00 
25.02 
105.00 


ROLL 
.  00 


FOLLOWING  OPTIONS  ARE  AVAILABLE 

ENTER  D  FOR  SEGMENT  POSITION  AND  ORIENTATIONS 

E  FOR  ENTERING  NEW  POSITION  AND  ORIENTATIONS 
Q  FOR  QUITTING 
Q 
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DO  YOU  ACCEPT  THIS 


INITIAL  CONDITION 


9 


Y 


DATA  DECK  IS  READY 
DATA  FILE  NAME  IS 


FOR  JOB  SUBMISSION 
CVSINPT 


Stop  -  Program  terminated. 


C: \ATB\UCP> 
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C : \ATB\USERCONF>atbposn 

ENTER  INPUT  FILENAME,  EXTENSION  . ain  IS  ASSUMED: 


driver 

ENTER  FILENAME  FOR  ALL  OUTPUT  FILES, 

EXTENSIONS  WILL  BE  ASSIGNED: 

driver 

ATBIV-4  COMPUTER  MODEL  -  CEA  VERSION 

ISTEP  =  0  TIME  =  0.00000  seconds 

ATB  -  INITIAL  POSITIONING  UTILITY 


SEGMENT  NUMBER  1  (  LT  )  LINEAR  POSITION  X=  -6.224  Y=  0.000  Z=  -4.92 


SEG 

YPRl 

YPR2 

YPR3 

LINEAR 

ACCELERATION  (G) 

ANG  ACCEL  (RAD/ S~ 2) 

X 

Y 

Z 

X 

Y 

Z 

LT 

0.00 

25.24 

0.00 

-0.810 

0.000 

-0.140 

0.0 

70.0 

0. 

CT 

0.00 

25.24 

0.00 

0.101 

0.000 

0.801 

0.0  - 

■141.0 

0. 

UT 

0.00 

25.24 

0.00 

-0.241 

0.000 

0.148 

0.0 

18.5 

0. 

N 

0.00 

12.73 

0.00 

-0.312 

0.000 

0.433 

0.0 

-61.2 

0. 

H 

0.00 

0.00 

0.00 

0.151 

0.000 

0.513 

0.0 

-5.7 

0. 

RUL 

0.00 

106.78 

0.00 

-0.501 

0.056 

-0.129 

-2.7 

-2.1 

4. 

RLL 

0 . 00 

50.68 

0 . 00 

0 . 503 

0.072 

-1.047 

2.4 

75.5 

4. 

RF 

0.00 

140.69 

0.00 

2.425 

-0.012 

-1.065 

-155.3  - 

171.8 

316. 

LUL 

0.00 

106.78 

0.00 

-0.501 

-0.056 

-0.129 

2.7 

-2.1 

-4. 

LLL 

0.00 

50.68 

0.00 

0.503 

-0.072 

-1.047 

-2.4 

75.5 

-4. 

LF 

0.00 

140.69 

0.00 

2.425 

0.012 

-1.065 

155.3  - 

171.8 

-316. 

RUA 

0.00 

25.02 

0.00 

-0.335 

0.000 

0.258 

0.0 

6.8 

-0. 

RLA 

0.00 

105.00 

0.00 

-0.071 

-0.002 

0.739 

0.2 

-28.6 

0. 

LUA 

0.00 

25.02 

0.00 

-0.335 

0.000 

0.258 

0.0 

6.8 

0. 

LLA 

0.00 

105.00 

0.00 

-0.071 

0.002 

0.739 

-0.2 

-28.6 

0. 

Enter 

X  ,  Y ,  Z  ,  L , 

E , J , ?  or 

SEG  &  YPR 

(i.e. 

"RUL2" ) 

to  change, 

0  to  END 

e 

(INERTIAL) 

( LOCAL ) 

U1  ARRAY  (INCH/ 

SEC**2 ) 

U2 

ARRAY  (RAD/ 

SEC**2 ) 

ANGULAR 

ACCELERATIONS 

SEGMENT 

X 

Y 

Z 

X 

Y 

Z 

1 

LT 

-309.49 

0.00 

-640.51 

0.00 

-32.26 

0.00 

2 

CT 

0.00 

0.00 

386.09 

0.00 

0.00 

0.00 

3 

UT 

19.90 

0.00 

376.80 

0.00 

-0.21 

0.00 

4 

N 

0.00 

0.00 

386.09 

0.00 

0.00 

0.00 

5 

H 

0.00 

0.00 

386.09 

0.00 

0.00 

0.00 

6 

RUL 

-196.38 

0.00 

-265.32 

0.00 

-396.09 

788.84 

7 

RLL 

0.00 

0.00 

386.09 

0.00 

0.00 

0.00 

8 

RF 

0.00 

0.00 

-3322.52 

-242.10 

-1475.97 

383.90 

9 

LUL 

-196.38 

0.00 

-265.32 

0.00 

-396.09 

-788.84 

10 

LLL 

0.00 

0.00 

386.09 

0.00 

0.00 

0.00 

11 

LF 

0.00 

0.00 

-3322.52 

242.10 

-1475.97 

-383.90 

12 

RUA 

0.00 

0.00 

386.09 

0.00 

0.00 

0.00 

13 

RLA 

0.00 

0.00 

386.09 

0.00 

0.00 

0.00 

14 

LUA 

0.00 

0.00 

386.09 

0.00 

0.00 

0.00 

15 

LLA 

0.00 

0.00 

386.09 

0.00 

0.00 

0.00 

Enter  X,Y 

,  Z  ,  L , E , J , ?  or 

SEG  &  YPR 

(i.e.  "RUL2" ) 

to  change, 

0  to  END  j 
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(INERTIAL) 

(INERTIAL) 

RELATIVE 

JOINT  FORCES  ( 

LBS) 

JOINT  TORQUES  (INCH 

LBS) 

ANGULAR  VEL 

JOINT 

IPIN 

X 

Y 

Z 

X 

Y 

Z 

(RAD/  SEC 

1 

P 

0 

-20.2 

0.0 

-66.5 

0.0 

0.0 

0.0 

0.0 

2 

w 

0 

-20.8 

0.0 

-65.4 

0.0 

0.0 

0.0 

0 . 0 

3 

NP 

0 

0 . 7 

0.0 

-6.0 

0.0 

0.0 

0.0 

0.0 

4 

HP 

0 

1 . 4 

0.0 

-4.6 

0.0 

0.0 

0.0 

0.0 

5 

RH 

0 

10 . 2 

2.0 

11.5 

0.0 

0.0 

0.0 

0.0 

6 

RK 

1 

10 . 1 

0.6 

-1.6 

-16.2 

0.0 

0.1 

0.0 

7 

RA 

0 

5.5 

0.0 

17.0 

0.0 

0.0 

0.0 

0.0 

8 

LH 

0 

10.2 

-2.0 

11.5 

0.0 

0.0 

0.0 

0.0 

9 

LK 

1 

10.1 

-0.6 

-1.6 

16.2 

0.0 

-0.1 

0 . 0 

10 

LA 

0 

5 . 5 

0.0 

17.0 

0.0 

0.0 

0.0 

0 . 0 

11 

RS 

0 

-1.9 

0.0 

-4.7 

0.0 

0.0 

0.0 

0 . 0 

12 

RE 

1 

-0.3 

0.0 

-1.2 

0.3 

0.0 

-0.2 

0.0 

13 

LS 

0 

-1.9 

0.0 

-4.7 

0.0 

0.0 

0.0 

0 . 0 

14 

LE 

1 

-0 . 3 

0.0 

-1.2 

-0.3 

0.0 

0.2 

0.0 

Enter  X, 

Y,Z,L 

,  E , J , ?  or 

SEG  &  YPR 

(i.e. 

"RUL2"  )  to 

change ,  0 

to  END 

1 

ATB  -  INITIAL  POSITIONING  UTILITY 


SEGMENT  NUMBER 

1  (  LT 

)  LINEAR  POSITION  X= 

-6.224 

Y=  0. 

000  z= 

-4.92 

SEG 

YPR1 

YPR2 

YPR3 

LINEAR 

ACCELERATION  (G) 

ANG  ACCEL  (RAD/S~2) 

X 

Y 

Z 

X 

Y 

Z 

LT 

0.00 

25.24 

0.00 

-0.810 

0.000 

-0.140 

0.0 

70.0 

0. 

CT 

0.00 

25.24 

0.00 

0.101 

0 . 000 

0.801 

0.0 

-141.0 

0. 

UT 

0 . 00 

25.24 

0.00 

-0.241 

0.000 

0.148 

0.0 

18.5 

0  . 

N 

0 . 00 

12.73 

0.00 

-0.312 

0 . 000 

0.433 

0.0 

-61.2 

0  . 

H 

0 . 00 

0.00 

0.00 

0.151 

0.000 

0.513 

0.0 

-5.7 

0. 

RUL 

0.00 

106.78 

0.00 

-0.501 

0 .056 

-0.129 

-2.7 

-2.1 

4. 

RLL 

0 . 00 

50.68 

0.00 

0.503 

0 . 072 

-1.047 

2.4 

75.5 

4  . 

RF 

0.00 

140 . 69 

0.00 

2.425 

-0 . 012 

-1.065 

-155.3 

-171.8 

316. 

LUL 

0 . 00 

106 . 78 

0.00 

-0.501 

-0.056 

-0.129 

2.7 

-2.1 

-4  . 

LLL 

0.00 

50.68 

0.00 

0.503 

-0 . 072 

-1.047 

-2.4 

75.5 

-4  . 

LF 

0.00 

140.69 

0.00 

2.425 

0.012 

-1.065 

155.3 

-171.8 

-316. 

RUA 

0.00 

25.02 

0.00 

-0.335 

0.000 

0.258 

0.0 

6.8 

-0. 

RLA 

0.00 

105.00 

0.00 

-0.071 

-0.002 

0.739 

0.2 

-28.6 

0. 

LUA 

0.00 

25.02 

0.00 

-0.335 

0 . 000 

0.258 

0.0 

6.8 

0. 

LLA 

0 . 00 

105.00 

0.00 

-0.071 

0.002 

0.739 

-0.2 

-28.6 

0. 

Enter 

X  ,  Y  ,  Z  ,  L  , 

E  ,  J ,  ?  or 

SEG  &  ’ 

yPR  (i.e. 

"RUL2 "  ) 

to  change 

,  0  to  END  rf 2 

Enter 

the  NEW 

VALUE  for 

(  RF 

),  YPR2  (Old  Value 

=  140. 

685)  141 
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ISTEP  =  0  TIME  =  0.00000  seconds 

ATB  -  INITIAL  POSITIONING  UTILITY 


SEGMENT  NUMBER  1  (  LT  )  LINEAR  POSITION  X=  -6.224  Y=  0.000  Z=  -4.92 

SEG  YPR1  YPR2  YPR3  LINEAR  ACCELERATION  (G)  ANG  ACCEL  (RAD/S~2) 


X 


LT 

0. 

.00 

25. 

.24 

0. 

.00 

-0 

.811 

CT 

0. 

.00 

25, 

,24 

0. 

.00 

0 

.101 

UT 

0, 

.00 

25, 

.24 

0. 

.  00 

-0 

.  241 

N 

0, 

.00 

12. 

.73 

0. 

.00 

-0 

,312 

H 

0. 

.00 

0, 

.00 

0, 

.00 

0 

.151 

RUL 

0. 

.00 

106. 

.  78 

0. 

.00 

-0 

.  504 

RLL 

0, 

.00 

50, 

.  68 

0  , 

.00 

0 

.  508 

RF 

0. 

.00 

141. 

,  00 

0. 

.00 

2 

.444 

LUL 

0. 

.00 

106, 

.  78 

0  , 

.00 

-0 

.  502 

LLL 

0. 

.00 

50. 

.  68 

0. 

,00 

0 

.503 

LF 

0. 

.00 

140, 

.  69 

0. 

.00 

2 

.425 

RUA 

0, 

.00 

25. 

.02 

0. 

.  00 

-0 

.  335 

RLA 

0  , 

.00 

105. 

,  00 

0, 

.00 

-0 

.071 

LUA 

0. 

.00 

25. 

.02 

0. 

,00 

-0 

.335 

LLA 

0, 

.00 

105, 

.00 

0, 

.00 

-0 

.071 

Y 

z 

X 

Y 

Z 

-0  , 

.001 

-0, 

.139 

0, 

.  1 

69, 

.  9 

0. 

0. 

.000 

0. 

.802 

0. 

.1 

-141. 

.1 

0. 

0  . 

.000 

0, 

.148 

0, 

.  0 

18, 

.  5 

0. 

0, 

.000 

0, 

.434 

0. 

,0 

-61, 

,1 

0. 

0  , 

.000 

0, 

.  513 

0, 

.0 

-5. 

.  7 

0. 

0. 

,058 

-0. 

.130 

-2. 

.8 

-1. 

,  9 

5. 

0  . 

.071 

-1, 

.063 

2. 

.7 

76, 

.2 

5. 

-0. 

.010 

-1. 

.115 

-157, 

.8 

-169. 

.  6 

319. 

-0, 

.057 

-0, 

,  129 

2. 

.7 

-2, 

.1 

-4. 

-0. 

.072 

-1. 

.047 

-2, 

.  4 

75. 

.  5 

-4. 

0. 

.012 

-1. 

.065 

155, 

.3 

-171, 

.8 

-316. 

0  . 

.000 

0. 

.258 

0, 

.0 

6. 

.8 

-0. 

-0. 

.002 

0, 

.739 

0. 

.2 

-28, 

.  6 

0. 

0, 

.000 

0. 

,  258 

0. 

.0 

6. 

.8 

0. 

0, 

.002 

0. 

.739 

-0. 

.2 

-28, 

,  6 

0. 

Enter  X , Y , Z , L , E , J , ?  or  SEG  &  YPR  (i.e.  "RUL2" )  to  change,  0  to  END  rf2 
Enter  the  NEW  VALUE  for  (  RF  ) ,  YPR2  (Old  Value  =  141.000)  140 

ISTEP  =  0  TIME  =  0.00000  seconds 


ATB  -  INITIAL  POSITIONING  UTILITY 

SEGMENT  NUMBER  1  (  LT  )  LINEAR  POSITION  X=  -6.224  Y=  0.000  Z=  -4.92 


SEG  YPR1  YPR2  YPR3  LINEAR  ACCELERATION  (G)  ANG  ACCEL  (RAD/S~2) 


X 


LT 

0  , 

.00 

25. 

.24 

0, 

.00 

-0. 

.809 

CT 

0. 

.00 

25. 

.24 

0  , 

.  00 

0  . 

,101 

UT 

0, 

.00 

25. 

.24 

0. 

.00 

-0. 

,241 

N 

0. 

.00 

12. 

.73 

0. 

,  00 

-0. 

,  312 

H 

0. 

.00 

0. 

.00 

0. 

.  00 

0  . 

.151 

RUL 

0  . 

.00 

106, 

.  78 

0, 

.00 

-0. 

.495 

RLL 

0, 

.00 

50. 

.  68 

0. 

.00 

0, 

.490 

RF 

0. 

.00 

140, 

.00 

0, 

.00 

2. 

,378 

LUL 

0, 

.00 

106, 

.78 

0. 

.00 

-0. 

.500 

LLL 

0. 

.00 

50. 

.  68 

0. 

.00 

0. 

.  503 

LF 

0. 

.00 

140, 

.  69 

0  , 

.00 

2, 

.425 

RUA 

0. 

.00 

25. 

,02 

0. 

.00 

-0. 

,335 

RLA 

0, 

,00 

105, 

.00 

0  , 

.00 

-0, 

.071 

LUA 

0. 

.00 

25, 

.02 

0, 

.00 

-0  , 

.335 

LLA 

0, 

,00 

105, 

.00 

0, 

.00 

-0. 

.071 

Y 

Z 

X 

Y 

Z 

0 . 002 

-0.141 

-0.2 

70.2 

-0  . 

0.000 

0.800 

-0.2 

-140.9 

0. 

0.000 

0.147 

0.0 

18.5 

0. 

0.000 

0.433 

0.0 

-61.2 

0. 

0.000 

0.513 

0.0 

-5.7 

0. 

0.053 

-0.126 

-2.4 

-2.7 

3. 

0.074 

-1.008 

1.8 

73.7 

4. 

-0.015 

-0.955 

-149.6 

-175.9 

310. 

-0.056 

-0.129 

2.8 

-2.1 

-4. 

-0.073 

-1.047 

-2.5 

75.4 

-5. 

0.012 

-1.065 

155.3 

-171.8 

-316. 

0.000 

0.258 

0.0 

6.8 

-0. 

-0.002 

0.739 

0.2 

-28.6 

0. 

0.000 

0.258 

0.0 

6.8 

0. 

0.002 

0.739 

-0.2 

-28.6 

0. 

Enter  X , Y , Z , L , E , J , ?  or  SEG  &  YPR  (i.e.  "RUL2"  )  to  change,  0  to  END  x 
Enter  the  NEW  VALUE  for  LINEAR  POSITION  -  X  (Old  Value  =  -6.224)  -7 
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ISTEP 


0 


TIME 


0.00000  seconds 


SEGMENT  NUMBER  1  (  LT 

SEG  YPR1  YPR2 


ATB  -  INITIAL  POSITIONING  UTILITY 

X=  -7.000 


LINEAR  POSITION 


Y  = 


0.000  Z=  -4.92 


YPR3 


LT 

0  . 

.00 

25 

.  24 

0. 

00 

CT 

0. 

,00 

25 

.24 

0. 

00 

UT 

0. 

.00 

25 

.  24 

0. 

00 

N 

0. 

,00 

12 

.73 

0. 

00 

H 

0  , 

.00 

0 

.  00 

0. 

00 

RUL 

0  . 

.00 

106 

.  78 

0. 

00 

RLL 

0, 

.00 

50 

.  68 

0. 

00 

RF 

0, 

,00 

140 

.  00 

0. 

00 

LUL 

0 

.00 

106 

.  78 

0  . 

00 

LLL 

0, 

.00 

50 

.68 

0. 

00 

LF 

0 

.00 

140 

.  69 

0. 

00 

RUA 

0 

.00 

25 

.02 

0. 

00 

RLA 

0 

.00 

105 

.00 

0. 

00 

LUA 

0 

.00 

25 

.02 

0. 

00 

LLA 

0 

.  00 

105 

.  00 

0. 

00 

2nter 

x,y,: 

Z,L, 

E,J, 

?  or 

SEG 

& 

LINEAR 

ACCELERATION  (G) 

ANG  ACCEL  (RAD/ S'- 2) 

X 

Y 

Z 

X 

Y 

Z 

-0.474 

0 . 002 

-0.035 

-0.2 

47.4 

-0  . 

0 . 523 

0.000 

0.765 

-0.2 

-141 . 6 

0  . 

0.095 

0.000 

0.143 

0.0 

23.5 

0  . 

-0.153 

0.000 

0.433 

0.0 

-38.8 

0  . 

0.154 

0.000 

0.489 

0.0 

-6.3 

0. 

-0 . 309 

0.144 

-0.081 

-6.7 

5.8 

18. 

0.546 

0.066 

-1.068 

11.6 

68.4 

15  . 

2.327 

-0.080 

-0.928 

-161.5 

-177.6 

282. 

-0.315 

-0 . 146 

-0.084 

7.1 

6.4 

-19. 

0.559 

-0.065 

-1.107 

-12.3 

70.1 

-16  . 

2.375 

0.077 

-1.038 

167.2 

-173.5 

-288. 

-0.189 

-0.002 

0.340 

0.1 

-3.6 

0. 

-0.084 

0 . 000 

0.787 

-0.2 

-24.3 

0. 

-0 . 189 

0.002 

0.341 

-0.1 

-3.6 

-0  . 

-0.084 

0 . 000 

0.787 

0.2 

-24.3 

-0. 

YPR  (i.e. 

"RUL2" ) 

to  change, 

o  to  : 

END  0 

C : \ATB\USERCONF> 
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Using  ATB  in  Collision  Reconstruction 

Wesley  D.  Grimes 

Collision  Engineering  Associates,  Inc. 


ABSTRACT 

The  Articulated  Total  Body  (ATB)  computer  program, 
sometimes  referred  to  as  the  Crash  Victim  Simulator  (CVS),  is  a 
powerful  tool  to  aid  in  studying  occupant  kinematics  in  motor 
vehicle  collisions.  There  are  many  options  available  within  the 
ATB/CVS  model  and  associated  utility  programs,  such  as 
GEBOD,  which  allow  an  analyst  to  model  specific  collisions  and 
occupants.  This  paper  discusses  ATB/CVS  as  a  tool  for  use  in 
collision  reconstruction.  Specific  examples  are  presented  in 
developing  a  crash  pulse  from  vehicle  simulation  programs  such  as 
EDSMAC,  SMAC,  HVOSM,  etc.  Techniques  are  also  presented 
for  modelling  other  moveable  objects  within  the  occupant 
environment,  such  as  a  seat  back,  steering  column,  or  intrusion  into 
the  occupant  compartment.  A  series  of  programs  to  aid  in  setting 
up  an  ATB  data  file,  the  CAL-3D  Users  Convenience  Package,  is 
also  discussed. 

INTRODUCTION 

The  Articulated  Total  Body  (ATB)[1]*  computer  program, 
also  known  as  the  Crash  Victim  Simulator  (CVS),  is  a  powerful 
simulation  tool  that  can  be  used  to  study  occupant  kinematics  in 
motor  vehicle  collisions.  However,  the  ATB/CVS  program  is 
generic  in  nature;  the  simulation  model  can  be  used  to  model 
occupants  in  motor  vehicle  collision,  occupants  being  ejected  from 
military  aircraft,  and  vehicle  rollover  kinematics,  to  name  a  few. 
This  generic  nature  may  complicate  the  setup  of  the  model.  There 
are  so  many  options  available,  a  casual  user  may  not  understand  the 
options  well  enough  to  make  a  proper  selection. 

There  are  also  many  options  that  were  introduced  in  earlier 
versions  of  the  program  that  have  been  superseded,  such  as  the  seat 
belt  being  replaced  by  the  harness  belt  However,  one  of  the  goals 
in  developing  the  ATB/CVS  model  has  always  been  to  maintain 
backward  compatibility  with  previous  data  files.  This  then  requires 
that  the  options  which  have  been  superseded  be  left  in  the  code  to 
ensure  that  older  data  files  will  still  produce  the  same  output 


*  Numbers  in  brackets  designate  references  found  at  the  end  of 
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Cards  A  -  Date  and  run  description,  units  of  input  and  output, 
control  of  restart,  integrator  and  optional  output 

Cards  B  -  Physical  characteristics  of  the  segments  and  joints. 

Cards  C  *  Description  of  the  vehicle  motion. 

Cards  D  -  Contact  planes,  belts,  air  bags,  contact 

(hyper)clIipsoids,  constraints,  symmetry  options, 
spring  dampers,  and  prescribed  forces  and  torques. 

Cards  E  -  Functions  defining  force-deflections,  stress-strains, 
inertial  spike,  energy  absorption  factor,  and  friction 
coefficients. 

Cards  F  -  Allowed  contacts  among  segments,  planes,  belts, 
airbags,  contact  (hyper)ellipsoids  and  harnesses. 

Cards  G  -  Initial  orientations  and  velocities  of  the  segments. 

Cards  H  -  Control  of  output  of  time  history  of  selected  segment 
motions,  joint  parameters,  wind  forces,  joint  forces 
and  torques,  total  body  properties,  and  injury  criteria. 

Cards  I  -  Control  information  for  plotter  output 


Figure  1  -  Outline  of  ATB  input 


There  is  a  subset  of  all  these  options  which  can  be  used  to 
model  occupants  in  automobile  collisions  without  the  burdensome 
task  of  understanding  ALL  the  other  options.  The  basic  options 
used  in  defining  the  occupant  as  well  as  the  vehicle  motion  will  be 
discussed  in  this  paperBvery  ATB  simulation,  including  those  used 
to  study  occupant  kinematics  in  a  motor  vehicle  collisions,  contains 
nine  (9)  groups  of  data,  identified  by  card  letters  A  through  I,  as 
shown  in  Figure  1.  The  input  file  describes  the  occupant,  the 
motion  of  the  vehicle,  and  the  possible  interactions  between  the 
occupant  and  the  environment  (typically  the  vehicle  interior). 
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OCCUPANT  DESCRIPTION 


One  of  the  tasks  in  setting  up  an  ATB  model  is  to  define  the 
parameters  for  the  occupant  The  occupant  is  typically  divided  into 
15  or  17  body  segments,  as  shown  in  Figure  2.  The  contact  surface 
for  each  segment  is  modelled  by  an  ellipsoid.  The  ellipsoid 
dimensions  along  with  inertial  properties  are  required  for  every 
segment  A  computer  program  has  been  developed  to  aid  a 
researcher  in  defining  these  parameters.  The  program  is  called 
GEBOD  (GEnerator  of BOdy  Data)[2].  Version  GEBODffl,  was 
released  in  1991,  with  version  GEBODIH.  1  as  an  update  in  April 
of  1993. 


The  GEBOD  program  can  be  used  to  produce  human  or 
anthropomorphic  test  device  (ATD)  models  for  use  in  the  ATB 
program.  There  are  three  different  ATD  types  which  can  be 
produced;  1)  sitting  Hybrid  III  (50  %tile),  2)  standing  Hybrid  HI 
(50  %tile),  and  3)  Hybrid  II  (50  %tile).  In  addition,  GEBOD  can 
read  a  file  with  specific  body  dimensions  supplied  by  the  user.  The 
program  then  uses  these  dimensions  to  generate  the  body  segment 
descriptions.  For  human  modelling,  the  program  prompts  the  user 
for  certain  specific  information  regarding  the  human  to  model, 
such  as  gender,  height,  weight,  etc.,  as  shown  in  Figure  3.  GEBOD 
then  uses  regression  equations  to  approximate  the  required 
properties  and  produce  a  file  in  the  format  required  by  ATB. 

Almost  any  size  occupant  can  be  modelled  using  GEBOD. 
Female  height  ranges  from  1 .446  to  1 .830  meters(56.93  to  72.05 
inches)  and  weight  ranges  from  378. 1  to  889.7  newtons  (85  to  200 
lbs).  There  are  similar  ranges  for  the  adult  male  occupant  model, 
with  the  height  ranging  from  1 .579  to  1 .972  meters  (62. 1 7  to  77.64 
inches)  and  the  weight  ranging  from  524.9  to  1 1 74.0  newtons  (118 
to  264  lbs).  A  child  model  is  also  available  which  permits  input  of 
age  from  2  to  20  years  old,  with  no  distinction  for  gender. 
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PROGRAM  GEBOD 

GEBOD  GENERATES  BODY  DESCRIPTION  DATA 
SUITABLE  FOR  INPUT  TO  THE  ATB  MODEL 
PLEASE  ENTER  A  DESCRIPTION  OF  THE  SUBJECT  (<60  CHARS.) 
Female  occupant,  62  inches  - 1 10  lbs 

ENTER  FILENAME  FOR  ALL  OUTPUT  FILES, 

EXTENSIONS  WILL  BE  ASSIGNED: 
sae-1 

1)  CHILD  (2 -20  YEARS) 

2)  ADULT  HUMAN  FEMALE 

3)  ADULT  HUMAN  MALE 

4)  USER-SUPPLIED  BODY  DIMENSIONS 
3)  SITTING  HYBRID  m  DUMMY  (50%) 

6)  STANDING  HYBRID  m  DUMMY  (50%) 

7)  HYBRID  n  DUMMY  (50%) 

ENTER  NUMBER  CORRESPONDING  TO  DESIRED  SUBJECT  TYPE 
2 

1)  WEIGHT 

2)  STANDING  HEIGHT 

3)  ALL  OF  THE  ABOVE 

ENTER  NUMBER  CORRESPONDING  TO 
PREDICTING  DIMENSIONS)  TO  BE  SUPPLIED 
3 

SELECT  UNITS  FOR  WEIGHT 

1)  LB. 

2)  N. 

3)  %-TTLE 


ENTER  VALUE  FOR  WEIGHT 
IN  THE  RANGE  85.00  ,  200.0 


SELECT  UNITS  FOR  STANDING  HEIGHT 

1)  IN. 

2)  M. 

3)  %-TILE 


ENTER  VALUE  FOR  STANDING  HEIGHT 
IN  THE  RANGE  56.93  ,72.05 


1)  FOREARM  AND  HAND  COMBINED 

2)  FOREARM  AND  HAND  SEPARATED 

ENTER  THE  NUMBER  CORRESPONDING  TO  THE  DESIRED 
LOWER  ARM  SEGMENTATION 
1 

SELECT  UNITS  FOR  OUTPUT 

1)  ENGLSH 

2)  METRIC 


Figure  3  -  GEBOD  prompts. 


Each  data  set  produced  by  GEBOD  must  be  checked  for 
reasonable  output  data.  The  program  uses  a  set  of  regression 
formulas,  and  sometimes  occupant  data  is  produced  that  is 
unreasonable  .  For  example.  Figure  4  shows  two  data  sets 
produced  for  a  small  female;  the  moments  of  inertia  for  the  center- 
torso  (CT)  are  negative.  Apparently,  the  regression  formulas 
produce  a  curve  fit  which  drops  below  zero  for  this  property  in  this 
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small  range  of  values.  This  problem  has  been  corrected  in  version 
GEBODIV,  which  should  be  distributed  in  early  1 995 .  Any  output 
used  from  previous  GEBOD  runs  modelling  a  small  female  should 
be  reviewed  carefully. 


Female  -  66 

Inches  in 

height  and  120  lbs 

CRASH  VICTIM 

PARAMETERS  (3-D) 

PRINCIPAL  MOMENT  OF  INERTIA 

SEGMENT 

WEIGHT 

(LB-SEC**2-IN ) 

I  SYM  PLOT 

(LB.  ) 

X  Y  Z 

1  LT  1 

14.351 

0.4332  0.2159  0.4541 

2  CT  2 

2.869 

0.0211  -0.0182  0.0041 

3  UT  3 

32.011 

1.9986  1.5555  1.1259 

Female  -  66 

inches  in 

height  and  110  lbs 

CRASH  VICTIM 

PARAMETERS  (3-D) 

PRINCIPAL  MOMENT  OF  INERTIA 

SEGMENT 

WEIGHT 

( LB-SEC**  2-IN ) 

I  SYM  PLOT 

(LB.  ) 

X  Y  Z 

1  LT  1 

11.915 

0.2882  0.0652  0.2191 

2  CT  2 

2.197 

-0.0138  -0.0480  -0.0568 

3  UT  3 

29.540 

1.7534  1.3643  0.9138 

Figure  4  -  Incorrect  properties  generated  by  GEBOD. 


VEHICLE  MOTION 

The  motion  of  the  motor  vehicle  can  be  described  using  one  of 
four  options,  as  shown  in  Figure  5.  For  studying  a  specific 
collision,  the  user  typically  uses:  Option  l  (half-sine  wave 
deceleration)  with  a  specified  speed  change;  or  Option  4  (spline 
fitted  position,  velocity,  or  deceleration)  using  results  of  a  vehicle 
simulation  program  such  as  EDSMAC[3],  SMAC[4],  or 
HVOSM[5].  The  method  described  herein  utilizes  the  results  of  an 
EDSMAC  simulation. 

Once  the  vehicle  simulation  has  been  completed,  the  resulting 
position  verses  time  data  is  used  as  input  for  the  ATB  program. 
The  velocity  or  acceleration  data  may  also  be  used,  it  simply 
depends  on  which  data  the  user  selects.  For  the  simplest 
description  of  motion,  the  position  and  orientation  data  are  the 
easiest  to  use. 

The  vehicle  position  data  which  is  output  from  most 
simulation  programs  is  in  feet.  This  data  must  be  converted  to 
inches  for  input  into  the  ATB  program.  The  EDSMAC  program 
stores  the  data  in  the  output  file  in  inches.  EDSMAC  then  converts 
the  data  to  feet  for  printout  A  proprietary  program  has  been 
Witten  which  will  read  the  data  in  the  EDSMAC  output  file  and 
generate  the  input  data  for  ATB.  This  utility  program  is  used  to 
quickly  produce  a  data  file  containing  the  EDSMAC  crash  pulse, 
in  ATB  format  This  data  file  is  then  inserted  into  the  ATB  input 
file  using  a  file  editor. 

The  spline  fit  routine  produces  a  data  look-up  table  of  vehicle 
accelerations,  either  directly  for  acceleration  entries  or  by 
erentiafing  position  and  velocity  data.  This  table  of  acceleration 
a.ta 15  used  by  the  program  just  as  if  the  user  had  entered  it 
usmg  Option  3  (six  degree-of-freedom  deceleration).  To  clarify, 


no  matter  what  type  of  data  is  input  (position,  velocity,  or 
acceleration),  the  spline  fit  routine  generates  a  table  of  acceleration 
data  that  is  then  used  by  the  program  integrator.  The  integration 
routines  in  the  ATB  model  will  then  integrate  the  acceleration  data 
in  these  tables  to  produce  the  velocity  and  position  of  the  vehicle. 
The  integration  routines  may  not  produce  exactly  the  same  position 
data  as  was  supplied,  depending  on  time  steps  of  integration,  type 
of  integration,  etc.  This  data  can  be  checked  using  printed  tables 
generated  with  the  H-cards  in  ATB. 

The  Rl,  R3,  R4,  andH.6  cards  should  be  used  to  verify  that 
the  motion  being  produced  by  the  ATB  model  is  what  the  user 
expects  and  that  the  motion  and  accelerations  of  the  vehicle  are 
reasonable.  The  H-cards  are  used  in  the  ATB  input  file  to  identify 
additional  time  history  tables  for  any  segments  or  joints.  These 
tables  may  include  linear  acceleration  (R 1),  linear  velocity  (H.2), 
and  linear  displacement  (H.3).  The  segment  angular  acceleration, 
velocity,  and  displacement  are  available  on  cards  H.4,  FL5,  and  H.6 
respectively.  Other  cards  (R7  through  RIO)  are  available  for 
joint  parameters,  wind  forces,  joint  forces  and  torques,  and  total 
body  properties.  The  ATB  program  will  also  calculate  the  Head 
Injury  Criteria  (HIC),  Head  Severity  Index  (HSI)  and  Chest 
Severity  Index  (CSI)  if  desired,  by  using  the  R 1 1  cards. 


Option  1 :  Half  sine  wave  deceleration  impulse  (NATAB  =  0) 

Option  2:  Tabular  unidirectional  deceleration  (NATAB  >  0) 

Option  3:  Six  degree  of  freedom  deceleration 
(NATAB  <  0  and  LTYPE  =  0) 

Option  4:  Spline  fit  position,  velocity  or  deceleration  data 
(NATAB  <0  and  LTYPE  >  0) 


Figure  5  -  Vehicle  motion  options  in  the  ATB  program. 

One  of  the  limitations  of  using  Option  4  for  specified  vehicle 
motion  is  that  only  101  data  points  may  be  entered.  In  order  to 
model  a  200  msec,  collision,  data  points  for  every  2  msec,  are 
output  from  EDSMAC.  For  a  long  ATB  run,  such  as  a  rollover, 
longer  time  intervals  may  be  used  or  other  vehicle  motion  options 
are  available  in  ATB. 

One  of  the  most  important  items  to  check  is  the  calculated 
vehicle  acceleration  produced  in  the  ATB  model.  Using  the  R 1 
and  R4  cards  the  acceleration  produced  in  the  ATB  model  can  be 
checked  against  the  test  data,  or  vehicle  simulation  data,  to  confirm 
that  the  approximate  acceleration  pulse  is  being  produced 
Similarly,  when  utilizing  acceleration  data  for  the  vehicle  motion 
data,  the  R3  and  R6  cards  should  be  used  to  confirm  that  the 
appropriate  displacements  are  being  produced  during  the  model 
integration. 

The  ATB  model  assigns  segment  numbers  to  every  object  that 
can  move  plus  one  for  the  ground  Thus,  with  15  occupant 
segments,  the  first  vehicle  is  segment  16  (NSEG+1).  Any  other 
vehicles  are  numbered  accordingly.  Air  bags  are  numbered  next, 
then  the  ground  is  numbered  as  the  highest  segment  number.  Thus, 
the  ground  (NGRND)  is  equal  to  the  number  of  occupant  body 
segments  (NSEG)  plus  the  number  of  airbags  (NBAG)  phis  the 
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number  of  vehicles  (NVEH)  plus  1 .  For  the  simplest  case  then, 
NSEG  is  15  for  the  occupant  body  segments,  NVEH  is  16  and 
NGRND  is  17.  Note  that  NGRND  can  usually  be  referred  to  as 
segment  number  zero  also. 

The  following  section  describes  suggested  values  on  the  ATB 
C -cards  when  using  a  crash  pulse  generated  by  EDSMAC.  The 
reader  should  have  the  ATB  User's  Guide  [1]  available  to  fully 
understand  the  variables  and  suggested  input  values.  In  order  to 
use  the  position  data,  generated  by  EDSMAC,  the  C -cards  should 
contain  the  following  data: 

Card  C.  1  -  Appropriate  title  describing  motion  data. 

Card  C.2.a  -  All  values  should  be  blank,  except: 

NATAB  -  Negative  value  of  the  number  of  acceleration 
points  to  be  produced,  this  includes  one  for  the  starting 
time.  For  example  (-201)  will  produce  201  spline  fitted 
data  points  describing  the  acceleration  of  the  vehicle.  The 
maximum  number  of  data  points  produced  is  50 1 . 

ATO  -  First  time  point  of  spline  fitted  data,  typically  this 
would  be  0.0  to  produce  a  table  that  starts  at  0.  seconds. 

ATD  -  Fixed  time  interval  (in  seconds)  for  spline  fitted 
data.  Typically  this  would  be  0.001  or  0.002  to  produce 
a  table  up  to  ABS(NATAB)*ATD  seconds.  For 
example,  200*0.001  would  produce  a  table  of  spline 
fitted  data  every  1  msec,  up  to  0.200  seconds,  while 
1 00*0.002  would  produce  a  table  of  data  every  2  msec, 
up  to  0.200  seconds. 

II  -  Flag  to  indicate  whether  the  prescribed  motion  is 
relative  to  the  ground  (0  -  default),  another  segment 
(11=2),  or  is  relative  to  another  segment,  but  with  an 
inverse  vector  (11=1).  The  inverse  option  is  used  to 
specify  the  motion  of  a  reference  system  relative  to  the 
vehicle  (useful  when  using  accelerometer  data). 


13  -  The  reference  segment  for  the  prescribed  motion  when 
thell  flag  is  not  zcto.  When  using  position  data  from  a 
vehicle  simulation,  this  should  be  zero  to  indicate  that  the 
positions  are  relative  to  the  ground 

MSEG  -  The  segment  number  associated  with  the 
prescribed  motion.  The  primaiy  vehicle  must  be  the  last 
motion  specified,  and  MSEG  must  be  zero  for  the  primary 
vehicle.  MSEG  being  zero  specifies  that  this  is  the  last 
motion  data  to  be  input 

Card  C.2.b  -  All  variables  should  be  blank  or  zero  except* 

LT  YPE  -  The  type  of  data  being  input  For  position  data, 
LTYPE  should  be  1,  for  velocity  data  LTYPE=2,  and  for 
acceleration  data,  LTYPE=3. 

LFIT  -  The  degree  of  polynomial  for  spline  fitting  the  input 
data.  For  position  data,  LFIT  must  be  2  or  3.  Usually  a 
value  of  3  will  produce  acceptable  results,  but  the  second 
derivative  (acceleration)  should  be  checked  by  using  the 
H.  1  and  H.4  cards. 

NPTS  -  Number  of  data  points  supplied,  the  maximum 
number  of  data  points  is  101 . 


Card  C.5  -  The  time,  in  seconds,  for  each  data  point  along  with 
the  position,  in  inches,  (X,YZ)  and  the  orientation  data,  in  degrees, 
(yaw,  pitch,  roll). 

Figure  6  shows  a  partial  output  listing  from  the  EDSMAC 
model  of  RICSAC  case  6.  The  motion  of  the  Volkswagen  Rabbit, 
vehicle  number  2,  is  modelled  in  this  example. 

The  corresponding  ATB  input  listing  is  shown  in  Figure  7. 
Recall  that  the  output  from  the  EDSMAC  run  is  in  feet,  while  the 
input  required  for  the  ATB  model  is  in  inches.  The  input  listing 
shows  the  options  for  position  input  as  described  above. 


RICSAC  Case 

#6  Chevelle 

vs .  Rabbit 

Time 

X  #2 

Y  #2 

PSI  #2 

V  tot  #2 

U-vel  #2 

V-vel  #2 

PSI-dot  #2 

Sec 

ft 

ft 

deg 

ft/sec 

ft/sec 

ft/sec 

deg/sec 

0.000 

12.17 

1.00 

120.00 

31.49 

31.49 

0.00 

0.00 

0.002 

12.14 

1.05 

120.00 

31.48 

31.48 

0.00 

0.00 

0.003 

12.12 

1.08 

120.00 

31.48 

31.48 

0.00 

0.00 

0.005 

12.09 

1.14 

120.00 

31.47 

31.47 

0.00 

0.00 

0.007 

12.06 

1.19 

120.00 

31.46 

31.46 

0.00 

0.00 

0.009 

12.03 

1.25 

120.00 

31.45 

31.45 

-0.00 

0.01 

0.010 

12.01 

1.27 

120.00 

31.45 

31.45 

-0.00 

0.01 

0.012 

11.98 

1.33 

120.00 

31.44 

31.44 

-0.01 

0.01 

0.014 

11.95 

1.38 

120.00 

31.42 

31.42 

-0.01 

0.02 

0.016 

11.92 

1.44 

120.00 

31.40 

31.40 

-0.02 

0.03 

0.018 

11.88 

1.49 

120.00 

31.38 

31.38 

-0.02 

0.05 

0.020 

11.85 

1.54 

120.00 

31.36 

31.36 

-0.03 

0.08 

Figure  6  -  Partial  EDSMAC  output  listing. 
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ADDITIONAL  VEHICLES  -  In  order  to  include  the  second 
vehicle  from  this  EDSMAC  run,  a  second  "vehicle"  segment  must 
be  defined.  This  vehicle  will  be  referred  to  as  a  "secondary" 
vehicle,  with  the  distinction  being  that  the  occupant  of  interest  is  in 
the  "primary"  vehicle.  In  the  ATB  data  input  file,  there  are  many 
references  to  the  "primary"  vehicle  and  some  variables  will  default 
to  this  segment,  such  as  occupant  reference  systems  for  initial 
positioning. 

The  secondary  vehicle  does  not  need  all  the  B-card 
information,  as  long  as  no  other  inertial  segments  will  be  attached 
to  it  .  The  exact  motion  will  be  specified  on  C-cards.  By 
specifying  the  next  available  segment  number  (MSEG  is  set  to 


NSEG+1),  the  ATB  program  sets  up  an  additional  segment  This 
segment  has  no  inertial  properties  and  only  contact  panels  or 
ellipsoids  can  be  attached  to  it  In  this  manner,  the  second  vehicle 
from  the  EDSMAC  run  can  be  specified.  A  contact  panel  is 
attached  to  investigate  the  interaction  of  the  occupant,  in  the 
primary  vehicle,  with  the  front  hood,  side,  etc.  of  the  second 
vehicle.  Figure  8  shows  a  partial  listing  of  an  ATB  input  file 
using  this  technique. 

In  this  data  file,  the  secondary  vehicle  is  segment  number  1 6, 
the  primary  vehicle  is  segment  number  17,  and  the  ground  is 
segment  number  1 8.  The  C-card  section  now  contains  two  sets 
of  motion  data  The  first  set  of  C-cards  is  for  the  secondary  vehicle. 


16  15  Largest  adult  female  CARD  B.l 

LT  1  31.001.44621.15071.9651  5.299  8.314  4.866  0.353  0.000  2.094  1  CARD  B.2 

0.00  0.00  0.0 


0.01  0.01  0.01  0.00  0.00  0.00  0.10  0.10  0.10  0.00  0.00  0.00CARD  B. 6 

EDSMAC  SIMULATION  OUTPUT  FOR  VEH  1  IN  RICSAC  #6  (ATB  SECOND  VEHICLE)  CARD  Cl 
0 •  0.  0.  0.  0.  0.  0.  0.  -200  0.  0.001  0  016CARD  C2A 

1  3  101  CARD  C2B 

.0000  .0000  .0000  0.0  .0000  0.0  0.0 

.2000  59.2789  4.6024  0.0  11.0669  0.0  0.0 

EDSMAC  SIMULATION  OUTPUT  FOR  VEH  2  IN  RICSAC  #6  CARD  Cl 

.  °*  0.  0.  0.  0.  0.  0.  -200  0.  0.001  0  0  0CARD  C2A 

*  3  101  CARD  C2B 

.0000  146.0000  12.0000  0.0  120.0088  0.0  0.0 

.2000  131.8055  66.8012  0.0  135.5376  0.0  0.0 

6°00000000  CARD  D1 


.0.  0.  0. 

1  3  101 

.0000  146.0000 

.2000  131.8055 

6  0  o 


12.0000 
66.8012 
0  0 


0.0  120.0088 
0.0  135.5376 

0  0  0 


1.000  1.000  1.000 

•12.0000  -13.0000  -13.3260 

0.0000  21.0000  0.0000 


0.0000  90.0000  0.0000 

1  0  17  0.0 


CARD  G1A 
0  17  CARD  G2A 
321  17CARD  G3A 


321  17CARD  G3 
CARD  H1A 


gure  8  Partial  listing  of  ATB  input  file  showing  parameters  used  for  "dummy"  vehicle. 
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with  MSEG  equal  to  16  on  card  C.2.  The  second  set  of  C -cards 
contains  motion  data  for  the  primary  vehicle,  with  MSEG  equal  to 
zero  on  card  C.2.  Recall  that  MSEG  being  zero  signals  the  end  of 
the  C-card  section;  thus  the  C-cards  for  the  primary  vehicle  must 
be  the  last  set  of  C-cards. 

This  type  of  ATB  run  is  useful  in  studying  the  interaction 
between  an  occupant  and  the  intruding  structure  of  another  vehicle. 
Care  must  be  taken  when  modelling  this  type  of  collision  using 
ATB  with  motion  from  a  simulation  such  as  EDSMAC.  The 
simulation  program  is  modelling  a  deformed  vehicle  while  ATB 
does  not  model  the  vehicle  deformation  at  all.  Thus,  when 
examining  the  interaction  with  the  opposing  vehicle,  the  panels 
used  to  define  the  vehicle  must  be  correctly  defined  to  model  the 
deformation  which  the  simulation  is  predicting.  If  this  is  not  done, 
the  striking  vehicle  will  appear  to  intrude  much  further  into  the 
struck  vehicle  than  the  simulation  is  showing,  and  the  results  wall 
not  be  correct. 

In  these  situations,  a  dynamic  crush  pattern  must  be 
established.  Additional  vehicle  panels  are  then  used  in  the  ATB 
model  to  depict  the  intrusion.  A  "secondary"  vehicle  is  set  up  with 
the  deformed  panels  attached  to  it  The  motion  for  the  panels  is 
then  defined  by  specifying  the  motion  of  this  secondary  vehicle. 
Then  using  the  deformation  pattern  from  EDSMAC  (at  the 
specified  time  from  start  of  impact)  and  performing  an  iterative 
series  of  runs,  the  time  of  interaction  between  the  occupant  and  the 
intruding  vehicle  can  be  confirmed. 

OTHER  MOVEABLE  SEGMENTS  -  Using  a  similar 
technique  of  a  "secondary"  vehicle,  other  moveable  interior 
components  may  be  modelled.  This  technique  involves  setting  up 
a  "dummy"  segment  in  the  B-cards.  Unlike  the  technique  described 
above  with  the  second  vehicle  from  EDSMAC,  another  segment 
must  be  used.  Any  time  a  segment  with  inertial  properties  will  be 
attached,  a  "dummy"  segment  must  be  used.  This  specification  is 
required  due  to  the  chain  setup  of  the  ATB  model.  The  model 
expects  each  segment  identified  in  the  B-cards  to  be  connected  to 
one  of  the  previously  defined  segments.  For  example,  the  sixth 
segment  specified  by  the  B.2  cards  is  expected  to  be  connected  to 
one  of  the  first  five  segments  by  some  type  of  joint.  Thus,  if  a  seat 
back  mass  is  defined  as  the  sixteenth  segment,  then  the  ATB 
program  expects  it  to  be  connected  to  one  of  the  fifteen  segments 
defined  previously  in  the  input  data  There  is  a  "null"  joint  option, 
but  this  indicates  that  the  segment  is  not  joined  to  anything,  which 
is  not  the  case  for  the  seat  back  being  modelled 


LOCKED  MASS  FOR 

SEAT 

BACK  REFERENCE 

0.  0. 

0. 

0.  0. 

0. 

0. 

0. 

1  2 
.0000 

6 

,0000 

.0000 

0.0 

.0000 

2.0000 

,0000 

.0000 

0.0 

.0000 

4.0000 

,0000 

.0000 

0.0 

.0000 

6.0000 

,0000 

.0000 

0.0 

.0000 

8.0000 

0000 

.0000 

0.0 

.0000 

10.0000 

,0000 

.0000 

0.0 

.0000 

To  achieve  the  appropriate  setup,  a  "dummy"  segment  is 
defined  as  the  sixteenth  segment  with  a  null  joint,  then  the  seat 
bade,  steering  column,  or  other  moveable  segment  is  segment  1 7 
and  is  attached  to  this  dummy  segment  with  the  joint  required 
This  dummy  segment  requires  values  for  mass,  moments  of  inertia, 
etc.  on  the  B-cards.  This  segment  will  be  used  for  reference  only 
and  will  be  rigidly  attached  to  the  primary  vehicle,  therefore  the 
actual  values  used  for  inertia  properties  are  not  important.  For 
simplicity,  all  the  required  values  are  set  to  unity  (1.0).  Blank 
cards  are  also  required  for  the  G.2  and  G.3  cards  because  the 
dummy  segment  will  be  a  reference  segment 

The  dummy  segment  is  rigidly  attached  to  the  primary  vehicle 
by  specifying  motion  cards  that  lock  the  position  throughout  the 
simulation,  as  shown  in  Figure  9.  By  attaching  the  dummy  segment 
at  the  vehicle  c.g.,  the  reference  system  stays  the  same  as  the 
primary  vehicle.  Note  that  the  dummy  segment  could  be  attached 
at  any  location  and  at  any  rotation  as  required  The  location  and 
orientation  of  this  dummy  segment  is  an  important  consideration  in 
simplifying  subsequent  references  using  position  data  of  the  other 
"linked"  segments. 

For  the  seat  back  model,  a  mass  segment  is  defined  with  the 
appropriate  seat  back  properties,  and  a  contact  plane  is  then 
attached  to  this  seat  back  mass  segment  The  contact  may  be 
modelled  using  segment  to  segment  modelling;  however,  a 
rectangular  contact  plane  more  accurately  models  the  geometry  of 
the  seat  back  contact  surface.  The  seat  back  mass  segment  is 
attached  to  the  dummy  segment  with  a  hinge  joint  at  the  proper 
location  on  the  B.3  cards,  and  the  joint  torque  characteristics  can 
then  be  defined  either  on  the  B.4a  cards  or  on  the  E.7,  F.5  cards. 
The  ATB  graphics  output  for  a  seat  back  model  in  a  rear  impact  is 
shown  in  Figure  10. 

When  modelling  a  steering  column,  the  lower  column  is 
attached  to  the  dummy  segment  with  a  ball-and-socket  joint.  The 
upper  column  is  attached  to  the  lower  column  with  a  slip  joint  The 
steering  wheel  and  hub  usually  are  then  attached  to  the  upper 
steering  column  using  another  ball-and-socket  joint  The  joint 
sp ring  characteristics  are  defined  on  the  B.4,  B.5  cards  using  data 
from  testing  or  other  published  data. 

When  setting  up  the  ATB  program  to  model  vehicle  intrusion 
or  interior  panel  movement,  the  dummy  segment  is  not  necessary. 
In  this  case  a  second  reference  segment  is  defined,  just  like  the 
EDSMAC  example  above.  This  reference  segment  can  then  be 
moved  with  respect  to  the  primary  vehicle  by  using  position  versus 
time  data  from  the  collision  simulation  or  test  data.  This  movement 


CARD  Cl 

-5  0.  2 .  21816CARD  C2A 

CARD  C2B 

0.0  0.0 

0.0  0.0 

0.0  0.0 

0.0  0.0 

0.0  0.0 

0.0  0.0 


Figure  9  -  A.TR  C-cards  used  to  lock  "dummy"  vehicle  to  primary  vehicle. 
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is  defined  on  the  C.2  cards  as  described  earlier,  and  appropriate 
contact  panels  are  attached  to  the  second  reference  segment  It  is 
extremely  difficult  to  model  any  interior  intrusion  with  detailed 
contact  surfaces.  The  ATB  model  limits  the  detail  achievable 
because  of  limits  on  number  of  segments,  contact  panels,  etc. 


However ,  the  ATB  model  does  have  the  capability  to  model  this 
scenario  with  enough  detail  to  provide  insight  into  occupant- 
interior  interaction  which  will  aid  any  reconstructionist  in 
understanding  the  collision. 

USER  CONVENIENCE  PACKAGE 

The  most  difficult  task  in  running  the  ATB  model  is 
assembling  the  data  file.  The  model  requires  force-deflection  data 
for  every  contact  that  may  occur,  joint  flexural  spring  and  damping 


characteristics  for  every  joint,  and  detailed  geometry  of  the  model 
environment.  Assuming  that  this  data  is  available  or  may  be 
obtained  through  testing,  the  task  of  generating  an  input  file  for  the 
ATB  program  is  still  very  difficult  A  series  of  programs  were 
written  in  the  early  1 980's  to  aid  in  assembling  the  data  file.  The 
programs  were  called  the  "Cal-3D  User  Convenience  Package",  or 
UCP[6].  These  programs  have  subsequently  been  ported  down  to 
IBM-PC  format 


The  UCP  consists  of  several  programs  that  aid  in  developing 
and  maintaining  a  database  of  vehicle  interiors,  occupant  models, 
ATB  force  deflection  curves,  etc.  Using  the  UCP  then  allows  for 
easier  and  quicker  assembly  of  an  input  data  file  for  the  ATB 
program.  Each  database  consists  of  the  actual  data  fields, 

structured  in  the  exact  format  required  by  ATB.  As  the  user  selects 

the  appropriate  data  to  include,  that  section  is  read  from  the 
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database  file  and  written  to  the  ATB  input  file. 


In  addition,  the  UCP  programs  keep  track  of  the  function 
numbers  and  will  assign  them  as  required  in  the  contact  section  of 
the  input  file  without  the  user  being  burdened  with  this  task.  The 
UCP  programs  are  then  used  to  assemble  the  ATB  input  file  by 
prompting  the  user  for  certain  information  regarding  expected 
contacts,  etc.  Figure  1 1  shows  a  partial  listing  of  the  program 
session. 

The  UCP  also  contains  a  program  that  attempts  to  position  the 
occupant  in  initial  equilibrium.  This  portion  of  the  package  is 
extremely  useful  in  establishing  the  initial  position  of  an  occupant 
in  the  model,  but  unfortunately,  the  routine  is  not  always  successful . 
This  portion  of  the  package  still  needs  additional  work  to  become 
completely  dependable  in  positioning  the  occupant 

DISCUSSION 

The  ATB  model  is  a  valuable  research  tool  for  use  in 
understanding  occupant  kinematics  in  motor  vehicle  collisions. 
The  model  can  be  difficult  to  use  due  to  the  generic  nature  and 
wide  variety  of  options  available.  There  are  a  set  of  general 
options  that  can  be  used  to  model  occupants  as  part  of  a  collision 
reconstruction  analysis.  The  ATB  model  can  be  successfully  used 
with  collision  simulation  programs  such  as  EDSMAC,  SMAC, 
HVOSM,  etc.  In  this  manner,  the  ATB  model  is  modelling,  as 
closely  as  possible,  the  collision  pulse  generated  by  the  collision 
simulation  programs. 

The  GEBOD  program  is  a  utility  program  which  generates 
body  description  data  that  can  be  used  in  the  ATB  model.  The 
GEBODIII  version  of  the  program  generates  invalid  data  for  small 
females  and  should  be  used  with  caution.  An  updated  version  of 
the  program  should  be  available  in  early  1 995. 

The  UCP  aids  in  setting  up  an  ATB  data  file  and  is  also  useful 
in  maintaining  databases  of  interior  characteristic  parameters.  The 
routines  in  the  UCP  used  to  establish  occupant  equilibrium  are 
useful,  but  require  additional  research. 

The  utility  programs  discussed  in  this  paper  are  extremely 
helpful  in  setting  up  an  ATB  data  set,  but  like  all  computer 
programs,  the  data  produced  must  be  analyzed  by  a  qualified  user 
to  assure  it  is  reasonable  and  that  it  makes  sense. 


PROGRAM  AVAILABILITY 

The  ATB  program  and  the  GEBOD  program  are  available  for 
a  distribution  fee  from  Armstrong  Laboratory  at  Wright-Patterson 
Air  Force  Base. 

DISCLAIMER 

The  cases  presented  in  this  paper  are  hypothetical.  The  data 
has  been  used  for  example  purposes  only  and  should  not  be  used 
in  any  real-world  analysis. 
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function  selection 

PLANE  NO.  1  SEAT  CUSHION 

FORCE-* DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 

(ENTER  0  FOR  NO  CONTACTS,  999  TO  SEE  LIST)  :  999 
FUNC.  NO.  FUNC.  NAME 

1  SEAT  BOTTOM  FDF 

2  SEAT  BACK  FDF 

3  STIFF  SURFACES  FDF 

4  DASH  -  FEMUR  FDF 

5  FRONT  DASH- FEMUR  FDF 

6  DASH  -  TORSO  FDF 

7  HEAD  -  WINDS HLD  FDF 

8  COLUMN  -  BODY  FDF 

9  STEERING  COL  E/A  FDF 

10  SEGMENT-SEGMENT  FDF 

11  HARNESS  (CABLE)  FDF 

12  HARNESS  (LAP)  FDF 

13  HARNESS  (SHLDR)  FDF 

14  HARNESS  FDF 

FORCE- DEFLECTION  FUNCTION 
ENTER  FUNCTION  NUMBER 

(ENTER  0  FOR  NO  CONTACTS,  999  TO  SEE  LIST)  :  1 

FUNCTION  SELECTION 
PLANE  NO.  1  SEAT  CUSHION 
INERTIAL  SPIKE  FUNCTION 
ENTER  FUNCTION  NUMBER 

ENTER  0  FOR  NO  INERTIAL  SPIKE,  999  TO  SEE  LIST)  : 

FUNCTION  SELECTION 

PLANE  NO.  1  SEAT  CUSHION 

ENERGY  DISSIPATION  (R)  FUNCTION 

ENTER  FUNCTION  NUMBER 

(ENTER  999  TO  SEE  LIST)  :  999 

FUNC.  NO.  FUNC.  NAME 

1  SEAT  BOTTOM  1 

2  SEAT  BACK  1 

3  STIFF  SURFACES 

4  DASH  -  FEMUR  1 

5  DASH  -  TORSO  1 

6  STEERING  COLMN 

7  SEG.  -  SEG. 

ENERGY  DISSIPATION  (R)  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  :  1 

FUNCTION  SELECTION 
PLANE  NO.  1  SEAT  CUSHION 
PERMANENT  SET  (G)  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  :  999 

FUNC*  N0*  FUNC.  NAME 

1  SEAT  1 

2  STIFF  SURFACES 

3  DASH  -  FEMUR 

:  DASH  -  TORSO 

5  STEERING  COLMN 

set  (G)  function 
uo«L^tion  number 

(ENTER  999  TO  SEE  LIST)  :  1 


FUNCTION  SELECTION 
PLANE  NO.  1  SEAT  CUSHION 
FRICTION  COEFFICIENT  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  :  999 

FUNC.  NO.  FUNC.  NAME 

1  FRICTION  MU- 1.0 

2  FRICTION  MU  -  0.75 

3  FRICTION  MU-0.5 

4  HARNESS  1 

5  HARNESS  2 

6  CONSTANT,  F-0.0 

FRICTION  COEFFICIENT  FUNCTION 
ENTER  FUNCTION  NUMBER 
(ENTER  999  TO  SEE  LIST)  :  2 


FORCE  DEFLECTION  FUNCTIONS 
..  FUNCTION  SELECTION  BY  PLANE 
PLANE  FUNCTION  NO. 

1  SEAT  CUSHION  1 

2  SEAT  BACK  2 

3  FLOOR  BOARD  3 

4  TOE  BOARD  3 

5  WINDSHIELD  7 


DO  YOU  WANT  TO  MAKE  ANY  CHANGES  (YES  OR  NO) 


FRICTION  FUNCTIONS 
FUNCTION  SELECTION  BY  PLANE 
PLANE  FUNCTION  NO. 

1  SEAT  CUSHION  2 

2  SEAT  BACK  2 

3  FLOOR  BOARD  1 

4  TOE  BOARD  1 

5  WINDSHIELD  3 


DO  YOU  WANT  TO  MAKE  ANY  CHANGES  (YES  OR  NO) 
PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE  1 
SEAT  CUSHION  :  3 

ENTER  SEGMENT  NOS.  FOR  CONTACT  WITH  THIS  PLANE 
CONTACT  NO.  1 
INPUT  SEGMENT  NO.  :  1 

CONTACT  NO.  2 
INPUT  SEGMENT  NO.  :  6 

CONTACT  NO.  3 
INPUT  SEGMENT  NO.  :  9 


PLANE  CONTACT  SELECTION 

ENTER  NO.  OF  SEGMENTS  CONTACTING  PLANE  2 
SEAT  BACK  :  3 

ENTER  SEGMENT  NOS.  FOR  CONTACT  WITH  THIS  PLANE 

CONTACT  NO.  1 

INPUT  SEGMENT  NO.  :  1 

CONTACT  NO.  2 

INPUT  SEGMENT  NO.  :  2 

CONTACT  NO.  3 

INPUT  SEGMENT  NO.  :  3 


PLANE 

NO.  TITLE 

1  SEAT  CUSHION 

2  SEAT  BACK 

3  FLOOR  BOARD 

4  TOE  BOARD 

5  WINDSHIELD 


SEGMENT  CONTACTS 
NO.  SEGMENTS 


1-LT 

$ 

6-RUL  , 

9-LUL 

1-LT 

f 

2-CT  , 

3-UT 

fi-RF 

9 

11-LF  , 

0-RF 

5-H 

t 

11-LF  , 

DO  YOU  WISH  TO  MODIFY  ANY  CONTACTS  (YES  OR  NO)  ?  N 
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Development  of  a 
Simulation  Database  for  the 
Advanced  Dynamic 
Anthropomorphic  Manikin  (ADAM) 
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Systems  Research  Laboratories,  Inc. 

Louise  A.  Obergefell,  Ph.D. 

Lt.  David  Johnson 

Armstrong  Laboratory 

Gregory  Thompson 
Systems  Research  Laboratories,  Inc. 


Introduction 


♦  Advanced  Dynamic  Anthropomorphic 
Manikin  (ADAM)  is  a  human  surrogate 
used  for  high-speed  ejection  testing 

♦  Full-scale  testing  is  expensive  and  time- 
consuming 

♦  Want  to  supplement  testing  with  computer 
simulations 

»  Provides  additional  data  not  available  from  tests 
»  Use  as  design  tool  before  prototype  testing 
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♦  ADAM  designed  to  reproduce 
human  characteristics 

♦  ADARTs  design  constraints: 

»  biofidelity 

»  joint  motion  measurement 
»  durability 

♦  Constraints  compromises 
differences  between  humans  and 
ADAM 


Overview 


♦  Approach 

♦  Data  Requirements 

♦  Testing  &  Modeling  Methods 

♦  Verification 

♦  Recommendations 
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Approach 


♦  Measure  ADAM’s  physical  characteristics 

♦  Analyze  and  format  for  input  into 
simulation  program 

♦  Verify  that  data  was  input  correctly 

♦  Validate  data  set 


Required  Data 

♦  Segment  and  joint  configuration 

♦  Segment  inertial  properties 

♦  Ellipsoid  dimensions  and  offsets 
4  Joint  locations 

4  Joint  torque  characteristics 
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ADAM  Model  Configuration 


Inertial  Properties 

Required  for  each  segment: 

»  Mass 

»  Center  of  gravity  location 
»  Principal  moments  of  inertia 
»  Orientation  of  principal  axes 

STandard  Automated  Mass  Properties 
(STAMP)  Testing  Procedures 
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Segment  Coordinate  System  and 
Ellipsoid  Dimensions 


♦  Z  semiaxis  derived  from  joint  locations 

♦  X  and  Y  semiaxes  derived  from 
circumference  and  diameter 
measurements 

♦  Segment  coordinate  system  origin  located 
at  CG 

♦  All  segment  dimensions  given  with 
respect  to  segment  coordinate  system 


Joint  Locations 

♦  Measured  during  STAMP  procedure  with 
three-dimensional  digitizer 


*  Referenced  to  both  adjoining  segments’ 
coordinate  systems 
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♦  ADAM  neck  taken  from  standard  Hybrid  III 
data  set 

»  Two  ball-and-socket  joints 
»  Each  has  three  angular  degrees  of  freedom 

♦  All  other  joints  use  Euler  joint  model 

»  Rotations  and  torques  applied  about  precession, 
nutation,  and  spin  axes 

»  Different  torque-angle  characteristic  for  each  rotation 
»  Individual  joint  axis  rotations  may  be  locked 
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Verification 


*  Static  &  Dynamic  Simulations 

»  Ensure  proper  data  input 

»  Ensure  proper  data  handling  by  ATB  program 


Summary 


♦  Large  ADAM  measured 

♦  ATB  data  set  developed 

♦  Initial  verification  conducted 
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Recommendations 


♦  Validate  against  horizontal  sled  and 
ejection  test  results 

♦  Develop  small  ADAM  data  set 

♦  Use  to  supplement  future  AL  tests 
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Opportunities  for  Casualty  Reduction  in  Rollover  Crashes 

K.  H.  Digges,  A.  C.  Malliaris,  and  H.  J.  DeBlois 

DeBlois  Associates 

United  States 

Paper  No.  94-S5-0-11 

1994  ESV  Conference 

EXCERPT 

CHARACTERIZATION  OF  THE  ROLLOVER  ENVIRONMENT 

Earlier  studies  [Terhune  1988,  Malliaris  1991]  suggest  that 
the  most  common  type  of  rollover  involves  a  vehicle  slipping 
sideways  and  incurring  tripping  acceleration,  which  induces  the 
rollover.  Terhune  found  that  70%  of  the  cars,  and  91%  of  the  vans 
and  pickups  had  tripping  induced  rollovers.  Malliaris  reported 
that  precrash  events  which  can  induce  lateral  slide  are  predictors 
of  rollover  incidence,  and  that  vehicle  angular  motion  about  the 
roll  axis  is  recorded  in  more  than  95%  of  the  rollover  cases. 

Analysis  of  NASS  rollovers  shows  that  two  thirds  of  all 
vehicles  in  rollovers  undergo  less  than  one  complete  revolution 
(three  or  less  quarter  turns) .  This  group  of  rollovers  accounts 
for  46.6%  of  the  harm  to  occupants.  The  remaining  54.4%  of  the 
harm  occurs  in  vehicles  which  undergo  one  or  more  revolutions.  The 
mean  precrash  speed  for  rollovers  is  50  mph,  compared  to  28  mph  for 
all  other  crashes.  For  rollovers  in  which  a  fatality  is  involved, 
the  mean  speed  is  63.4  mph,  compared  with  45.3  mph  for  all  other 
fatal  crashes  [Malliaris  1991] . 

Estimates  of  the  roll  rate  were  made  by  analyzing  the 
trajectory  of  vehicles  from  140  cases  reported  in  NASS.  The  cases 
were  selected  from  the  900  rollovers  in  1989-90  NASS  on  the  basis 
of  extent  of  damage.  All  cases  in  the  file  with  severe  damage 
(Collision  Damage  Classification  of  4  or  more)  were  included  in  the 
study.  The  mean  precrash  speed  for  these  cases  was  56  mph,  and  the 
mean  roll  rate  was  between  1  and  2  rev/sec.  In  89%  of  the  cases, 
the  rollover  was  the  most  harmful  event.  The  tripping  force  was 
the  most  harmful  event  in  most  of  the  remaining  cases.  The  roll 
direction  was  43%  clockwise,  46%  counterclockwise,  and  11%  unknown. 

SIMULATION  OF  OCCUPANT  KINEMATICS  IN  ROLLOVERS 

Testing  of  countermeasures  in  rollover  has  been  extremely 
limited  to  date.  Laboratory  testing  is  limited  by  the  lack  of 
facilities  which  can  subject  a  dummy  to  the  rollover  crash 
environment.  There  is  no  rollover  equivalent  to  the  crash  test 
sled  which  is  used  for  developing  safety  systems  for  planar  impact 
environments.  A  great  deal  of  useful  information  has  been  gained 
from  rollover  testing  of  complete  vehicles.  During  the  past  ten 
years,  rollover  tests  conducted  by  the  Department  of  Transportation 
have  been  of  two  types. 

The  first  type  of  rollover  is  induced  by  a  median  barrier 
designed  to  redirect  vehicles  back  into  the  roadway.  However,  as 
the  angle  and  speed  of  engagement  with  the  barrier  are  increased, 
a  rollover  may  result.  The  roll  rate  observed  when  a  .midsize 
vehicle  traveling  at  60  mph  engages  the  median  barrier  is  about  1 
rev/ sec. 

The  second  type  of  test  is  a  staged  rollover  induced  by 
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ejecting  the  vehicle  from  a  moving  a  test  ij?  to -the  direction  of- 
contains  platform  which  is  hl"9ed  P?rP=n^“ifrat  a„  angle  (usually 

travel.  The  vehicle  is  placed  on  the  platform  at  an  ang  ^^ 

■£trj?z,sjrx ■“  STS  j,  “rj: 

generally  less  than  1  rev/ sec.  -ra-hes  permits  precise  and 

?°Elf  control  o?9  the  r'oU^r  envfron^  and  "permits  the 
repeatable  control  or  rne  ruj.xu  low  COst.  Computer 

study  of  a  wide  range  °f  an^occupants^in  rollover  crashes  have  been 
simulations  of  vehicles  and  ?°cupants  and  published  in  the 

developed  under  sponsorship^  _  .  These  papers  describe 

literature  [Obergefell  1986  Smith  iWl^These  pape^  ^ 

models  which  have  been  va  these  rollover  test  types, 

described  above.  However,  in  1 b?^  ro^  relatively 

the  tripping  acceleration  whfdd  ?d  this  ® /er  introduce  tripping 

scenarios  for  rollover  Simula  10  occupant  motion  during 

initial  goal  of  the  simulation  is  to  explore  occupant^  ^  ^  pre- 

the  initial  phase  of  the  .rolloY®.r  *  .  tripping  and  launching 

roll  conditions  including  ■  «*iad“g'  ‘“t  Unear  and  angular 
accelerations.  It  includes  t  d  first  impact  with  the 

motion  of  the  vehicle,  ending ’“rted  S  this  paper,  the  vehicle  is 

!ni t i al lyF°^l iding 1  sideway st  ^lefaM 

SS&TS  delta  V^.ItTs^liidtn6  c^n^ion  with  angular 
acceleration  over  a  period  of  120  ms. 

DISCUSSION  OF  EJECTION  CONTROL 

Increased  restraint  useis  a  read^Y  *V**^^^ 

reducing  ejections.  Jf.  e  provides  majo/ benefit  in 

occupant  containment  within  the  P  ,  ] .  •  871.  Based 

‘sr  .r 

windshield.  .  .  +._  avamine  the  side  window  ejection 

The  simulation  selected  earlier.  Variation 

opportunity  was  the  tripped ro  the  seVerity  of  the  tripping 

of  two  variables  were  mvestig  •  .  rate  The  occupant  was 

event  (specified  by  deltaV) ,  and  tte  P^  ^ 

motion  *wai^  SSETi'S 

ESS-? 

Ejections^re  not  produced  dur^\r°^OVe^wever,5  a^thi Gripping 
when  the  tripping  delta  was  only  ?;  “Ph‘  th  head  resulted  at  0.5 

severity  increased  to  5  rnph,  e^' ffL^iid  at  1  0  rev/ sec.  At  higher 
rev/sec,  and  complete  ejection  resulted  a  l-Ojev/  both  ;oll 

tripping  severities,  complete  ejection  was  prou 

rates . 
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.  -  ...  emulations  with  election  resistant  glazing  m 

the  sWeTindow  were  Bade  The 

i^nTHrHSS  sss^-sss 

glazing  by  the  dummy  were  well  below  the  strength  ox 

windshield  interlayers.  testing  of  ejection  resistant 

NHTSA  has  reported  successful  testing^ ^  evaluated 

side  windows  and  windshields  [Cla  alazing  by  laboratory  tests  and 
the  penetration  resistance  of  the  qMing  w  ^  vehicles.  The 

the  rollover  Perfo™“cfeorbyec1re|fiotneSresistant  glazings  included 
stated  design  goais  f°r  egection  and  maintaining  integrity 

resisting  a  40  lb  ball  impact  at  20  i mph,  and  m  vehlcles  with 

during  a  rollover.  Roll  windows  have  been  reported, 

experimental  ejection  resistant  side  ^indows^n^  ^  ^  tests.  In 

The  integrity  of  the  side  gla  ^  quarter  turns  ranged  from 

??oeif  and^the; vXcle^de^rBation  at  the  right  A-pillar  ranged 

fr°m  Additional"research  and a"aly^|^attd  "by^atva^cefgKzSg 
technology  "and  "to  'estimate*  the  "benefits  which  might  be  achieved. 


TABLE  III 

EJECTION  PATHS  IN  ROLLOVER  CRASHES 
an  UCHT  VEHICLE  OCCUPANTS 


Figure  2.  Rollover  Modal.  Head  Ejection 
0.6  Rev /Sec;  6  MPH  Tripping  Deita-V 
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The  Use  of  ATB/DYNAMAN  in  Injury  Biomechanics 


INTRODUCTION 

The  Articulated  Total  Body  (ATB)  program 
was  originally  designed  by  CALSPAN  with  the 
intent  of  improving  occupant  safety  and 
developing  modern  safety  standards  for  the 
automotive  industry.  In  recent  years,  the 
versatile  software  has  been  applied  to  the 
analyses  of  injury  biomechanics  resulting  from 
various  accidents. 

To  date,  ATB  has  been  used  to  model  the 
kinematics  and  dynamics  of  occupants  in  a 
much  wider  range  of  accident  scenarios, 
including  estimating  joint  forces  and  moments 
and  segment  accelerations  and  contact  forces. 
If  the  estimated  motions  and/or  loads  of  an  ATB 
analysis  are  compared  to  critical  injury- 
producing  values  which  have  been  estimated  or 
experimentally  obtained  from  the  literature,  then 
one  can  estimate  the  likelihood  of  various 
injuries  in  a  particular  accident.  In  some  cases, 
detailed  vehicle  motions  are  not  known  a  priori 
and  must  be  estimated  along  with  the  occupant 
kinematics  so  that  potential  injuries  may  be 
analyzed.  Three  examples  of  ATB  simulations 
are  discussed  which  demonstrate  how  ATB  can 
be  used  to  perform  an  injury  analysis  for 
different  vehicular  accidents. 

FRONTAL  COLLISION  EXAMPLE 

A  sports  car  underwent  a  35  mph 
longitudinal  Delta-V  in  a  head-on  collision  with 
the  side  of  a  swerving  van.  The  unbelted  front 
seat  passenger  of  the  automobile  sustained  a 
neck  fracture  resulting  in  quadriplegia,  a 
forehead  laceration,  and  loss  of  consciousness. 
The  ATB  program  was  used  to  investigate  the 
effectiveness  of  a  3-point  seat  belt  in  this 
scenario. 

An  ATB  model  of  the  car  interior  was 
created  based  on  measurements  of  an  exemplar 
vehicle.  Geometric  and  inertial  properties  of  the 


Tara  Khatua,  Louis  Cheng, 

Robert  Fijan,  Renate  Schmidt-Hargrave 

Failure  Analysis  Associates,  Inc. 


occupant  were  created  from  height,  weight,  and 
gender  data  automated  in  the  GEBOD  program. 
Vehicle  stiffness  data  were  taken  from  the  data 
set  supplied  along  with  the  ATB  program  [1-3]. 

The  simulation  of  the  unbelted  occupant 
produced  joint  forces  and  moments  consistent 
with  the  resulting  injuries,  while  a  simulation 
which  included  a  3-point  seat  belt  showed  that 
the  occupant  would  have  experienced  no  head 
or  knee  contact  with  the  vehicle  interior  and 
would  not  have  been  seriously  injured. 

FORKLIFT  ACCIDENT  EXAMPLE 

In  another  application,  ATB  was  used  to 
analyze  injuries  associated  with  a  forklift 
accident.  In  this  case,  ATB  was  used  to  not  only 
estimate  the  dynamics  of  the  occupant  of  the 
forklift  but  also  of  the  forklift  itself.  In  most 
automobile  accidents,  the  mass  of  the  vehicle  is 
at  least  an  order  of  magnitude  greater  than  the 
mass  of  the  occupant,  and  hence  one  can 
assume  that  the  vehicle  dynamics  are  not 
significantly  affected  by  the  occupant  kinematics. 
In  this  case,  however,  the  seat  was  attached  to 
the  vehicle  through  a  hinge,  and  hence  the 
seat's  motion  needed  to  be  considered. 
Therefore,  a  coupling  dynamic  analysis  including 
the  forklift  and  occupant  was  performed. 

The  geometry  of  the  forklift  chassis,  tires, 
seat,  and  the  rollover  protection  structure 
(ROPS)  in  the  ATB  model  were  obtained  by 
measurements  of  an  exemplar  forklift.  The 
center  of  gravity  and  inertia  data  were  computed 
using  the  geometric  data  and  standard  formulae. 
Seat  properties  were  obtained  by  performing 
compression  tests. 

Since  the  head  injury  parameters  were  the 
subject  of  concern,  care  was  taken  to  determine 
this  parameter  rigorously.  A  beamlike  structural 
stiffness  was  added  to  the  head  stiffness 
obtained  from  the  biomechanics  literature  [4-6]. 

A  dynamically-coupled  ATB  analysis  of  the 
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forklift  was  performed  to  analyze  forces  between 
the  subject's  head  and  a  ROPS  (rollover 
protective  system).  Results  showed  that  the 
occupant's  head  impacted  the  ROPS  structure 
before  the  ROPS  structure  impacted  the  ground, 
and  hence  the  injury  was  not  caused  by  forces 
to  the  head  from  a  deforming  ROPS. 

BUGGY  ROLLOVER  EXAMPLE 

A  final  example  involves  simulating  a  dune 
buggy  in  a  rollover  accident.  As  in  the  second 
example,  this  case  required  simulating  both  the 
occupant  and  vehicle  kinematics  using  ATB. 

The  ATB  model  was  created  based  on 
measurements  of  an  exemplar  vehicle,  and 
inertial  properties  were  determined  from 
experimental  tests.  Witness  statements  were 
used  to  estimate  the  starting  point  of  the  roll  and 
the  final  rest  position.  The  simulation  of  the 
roller  provided  understanding  of  the  inertial  loads 
imparted  to  the  vehicle  occupants  and  the  seat 
belt  restraint  system. 

SUMMARY 

By  using  ATB  to  calculate  joint  moments  or 
segment  contact  forces  and  combining  these 
estimates  with  critical  values  from  the  literature, 
an  injury  analysis  can  be  performed.  Three 
examples  have  been  discussed  which 
demonstrate  many  of  the  important  issues 
involved  in  performing  an  injury  analysis  .via 
computer  simulation  using  the  ATB  program. 
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The  Use  of  ATB/DYNAMAN  in 
Injury  Biomechanics 
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-  General  purpose  rigid-body  dynamics 
programs  are  very  tedious  and  time 
consuming  for  occupant  modeling. 


Introduction 
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this  program  can  still  be  used  as  a 
general  rigid-body  dynamics  (open  loop) 
program  with  impact  capabilities. 
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sand  dunes 


History  of  ATB  use  at  FaAA 
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Examples 
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+  Neck  fracture  resulting  in  quadriplegia 
+  Forehead  laceration  and  loss  of 
consciousness 


Effectiveness  of  seatbelt  usage 
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-  Vehicle  Interior  stiffness  data  :  Most  of  the 
data  were  taken  from  data  set  supplied 
along  with  the  ATB  program. 


Effectiveness  of  seatbelt  use 
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-  Unbelted  occupant  model :  The  model  for 
the  unbelted  occupant  is  shown  in  the  next 
slide. 
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Effectiveness  of  seatbelt  use 
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injury.  This  correlates  very  well  with  the 
passenger’s  loss  of  consciousness  and 
lack  of  permanent  head  injury. 


Effectiveness  of  seatbelt  use 
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slide. 


Effectiveness  of  Seatbelt  Use  -  Belted  Occupant 


Effectiveness  of  seatbelt  use 
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Forklift  Rollover  Simulation 


Forklift  rollover  simulation 
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after  it  bent  upon  hitting  the  ground. 

(2)  Head  was  injured  by  the  ROPS  before 
ROPS  struck  the  ground. 


Forklift  rollover  simulation 
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Forklift  rollover  simulation 
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Therefore,  a  combined  (dynamically 
coupled)  model  of  the  forklift  chassis,  seat, 
and  the  occupant  is  required. 


Forklift  rollover  simulation 
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Forklift  rollover  simulation 
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Forklift  rollover  simulation 
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Dune  Buggy  Rollover  Simulation 


Dunebuggy  rollover  simulation 
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Complexiety  of  the  modeling  :  Occupants 
are  interacting  with  dunebuggy  structural 
components  and  the  ground. 


Dunebuggy  rollover  simulation 
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Dunebuggy  rollover  simulation 
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Use  of  DYNAMAN  for  Design  and  Development  of  Anthropometric  Test 

Devices 

T.  Shams  and  Y.  Zhao,  GESAC,  Inc. 


1.  Introduction 

DYNAMAN  is  being  used  extensively  as  a  tool  to  aid  in  the 
design  of  a  number  of  components  for  an  enhanced  anthropometric 
test  device  which  is  currently  under  development.  The  ATD,  which 
is  called  the  Trauma  Assessment  Device  (TAD) ,  will  have  a  number 
of  enhanced  components,  with  improved  anthropometry  and 
biofidelity  [1].  These  include  the  thorax,  abdomen,  neck, 
head/face,  and  lower  leg  assembly.  In  this  presentation,  we  will 
provide  a  number  of  examples  where  DYNAMAN  has  been  used  to  model 
several  components  and  their  responses  under  impact  forces  or 
accelerations.  In  particular,  the  modelling  of  a  new  abdomen 
assembly  and  a  new  neck  will  be  described. 

The  DYNAMAN  models  are  being  used  to  provide  a  fast, 
practical  procedure  for  selecting  the  basic  characteristics  of 
different  hardware  parts  of  a  component,  such  its  geometry  and 
stiffness.  The  models  allow  one  to  vary  the  important  design 
parameters  to  under  the  sensitivity  of  the  response  to  these 
parameters.  In  some  cases,  the  modelling  is  being  used  as  a 
preliminary  procedure  for  identifying  these  characteristics 
before  proceeding  with  more  advanced  modelling  using  finite 
elements.  But  in  our  experience,  we  have  found  that  for 
materials  such  as  rubbers  and  foams,  which  are  somewhat  difficult 
to  characterize  within  the  context  of  a  finite  element  model  is 
more  easily  defined  using  the  lumped  parameters  used  in  DYNAMAN, 
where  only  the  average  response  is  required. 

The  first  step  in  all  the  modelling  exercises,  is  to  select 
the  modelling  elements  (e.g.  ellipsoids,  joints,  spring-dampers, 
airbags,  belts),  that  will  best  describe  a  particular  portion  of 
a  dummy  component.  The  appropriate  geometry  associated  with  a 
modelling  element  is  next  determined.  The  geometry  is  usually 
guided  by  the  actual  physical  features  of  the  part,  though  in 
some  cases,  compromises  will  need  to  be  made  in  order  to  fit  a 
lumped  parameter  to  what  is  really  an  extended  object.  The  next 
objective  in  the  modelling  is  to  select  the  response  parameters, 
which  would  make  the  component  behave  in  a  way  similar  to  what  is 
suggested  by  the  design  goals.  These  may  include  the  force- 
deflection  characteristics  of  an  impact  surface,  elastic  and 
viscous  properties  of  a  spring  damper,  etc. 


2.  Modelling  of  the  Abdomen 

The  current  abdomen  is  a  frangible  insert  made  of . styrofoam, 
which  was  originally  designed  to  be  used  with  the  Hybrid  III 
dummy  [2].  It  has  "human-like"  force-deflection  characteristics 
when  loaded  by  belts,  but  can  only  measure  the  maximum 
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deformation  and  has  to  be  replaced  after  each  test.  A  number  of 
researchers  have  examined  the  possibility  of  using  a  foam-filled 
bag  as  the  surrogate  structure  for  the  abdomen  [3,4]. 

As  part  of  a  detailed  analysis  of  the  possible  alternatives 
that  could  be  used  for  the  design  of  a  new  abdomen,  an  lnitia 
model  was  developed  with  DYNAMAN  which  described  a  structure 
which  included  both  the  effects  of  the  foam  and  of  the  enclosed 
air.  The  enclosed  air  was  modeled  because  one  of  the  possib 
choices  for  the  foam  was  to  selected  an  open-cell  foam  where  a 
significant  percentage  of. the  volume  of  the  foam  is  taken  by  air. 
During  the  compression  of  a  such  a  structure,  resistance  will  be 
offered  by  both  the  elastic/ damping  characteristics  of  the 
material  of  the  foam,  and  by  the  "air  spring"  effect  of  the 
compressed  volume  of  air. 

The  initial  steps  in  the  development  of  the  model  consisted 
of  the  following: 

1.  Determination  of  the  approximate  geometry  of  the  structure. 
This  included  the  dimensions  available  for  assembly  in 
height,  weight,  and  breadth.  The  dimensions  were  dictated 
partly  by  the  space  available  in  the  lower  section  of  th_ 
thorax  between  the  lower  ribcage  and  the  pelvis.  The  main 
design  restrictions  were  the  area  of  the  abdomen  that  would 
be  exposed  in  the  front.  This  should  be  reasonably  simil 

to  the  exposed  area  of  a  human  abdomen  which  may  be  impacted 

2.  Determination^of^the  constraint  conditions  on  the  structure. 
Alternatives  included  structures  which  could  rotate  about 
the  lumbar  spine  (in  the  pitch  direction)  and  those  where 
the  abdomen  was  more  rigidly  attached  between  the  spine  and 

the  front  of  the  lower  thorax.  #  ^ 

3.  Determination  of  the  different  impact  situations  that  would 
need  to  be  modeled.  At  the  preliminary  stage  our  main 
interest  was  to  model  the  effects  of  a  rigid  rod  impacting 
the  upper  portion  of  the  abdomen  (similar  to  a  steering 
wheel  impact)  and  that  of  a  lap  belt  impacting  the  lower 

abdomen . 

The  basic  components  of  the  model  were: 

1.  A  spring-damper  system  which  modelled  the  foam 

characteristics  within  the  bag.  .  .  - 

2.  An  airbag  which  modelled  the  air  spring  characteristics  of 

3  AhSmall'mass  segment  (with  associated  ellipsoid)  which 

modened  the  initial  and  stiffness  characteristics  of  the 
bag  cover  which  would  be  in  immediate  contact  with  t 
impacting  object. 
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Figure  1.  General  layout  of  abdominal  segment  within 
the  dummy. 


Figure  2 .  Initial  setup  of  spring-damper  plus  airbag  model 
of  abdomen. 

Figure  1  shows  the  layout  of  the  area  of  the  dummy  where  the 
abdominal  segment  will  reside.  The  principal  attachment  is  with 
the  spinal  column  and  the  figure  shows  that  the  upper  abdominal 
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area  overlaps  the  lower  ribcage  and  the  lower  abdominal  area 
which  sits  on  top  of  the  pelvis.  Figure  2  shows  the  basic 
abdominal  model  which  consists  of  a  spring-damper,  and  airbag,  a 
mass  segment  (the  front  ellipsoid)  which  represents  the  cover. 

The  spring-damper  attaches  the  lower  torso  segment  (lower  spine) 
with  the  abdomen  cover  segment.  An  ellipsoidal  airbag  was 
defined  (though  a  cylindrical  bag  could  also  be  defined)  and  wa 
attached  to  the  lower  torso  segment.  The  airbag  was  inflated  in 
such  a  way  that  it  was  fully  inflated  (at  atmospheric  pressure) 
before  any  impact  could  occur  (DYNAMAN  cannot  generate  a 
preinflated  bag,  but  has  to  use  a  mass  flow  rate  function  to  feed 
gas  into  the  bag.)  The  initial  model  analyzed  the  impact  of  a 
steering  wheel  (represented  by  an  ellipsoid)  at  an  angle  of  45 
with  the  bag  cover  segment. 

NOTE:  Some  care  has  to  be  taken  in  defining  the  contact 

region  of  the  bag  with  the  impactor,  since  the 
penetration  algorithm  used  can  fail,  once  the  center  of 
the  penetrating  segment  passes  the  center  of  the  ag. 

Once  the  basic  model  was  operational,  the  whole  bag  assembly 
was  then  inserted  within  the  Hybrid  III  dummy  (replacing  the 
existing  lower  torso) .  This  model  was  exercised  using  both 
impacts  with  a  steering  wheel  and  with  a  lap  belt.  Figure  3 
shows  the  setup  of  the  whole  dummy  with  the  abdomen  component 
within  it,  and  with  a  lap  belt  attached.  Figure  4  shows  the 
picture  of  the  dummy  at  the  time  of  maximum  penetration  of  the 
belt  into  the  abdominal  bag. 

The  objective  of  the  model  was  to  estimate  the  possible 
ranqe  of  values  for  the  foam  stiffness  and  of  the  airbag  volume 
which  would  result  in  responses  similar  to  that  seen  for  humans. 
Unfortunately  there  is  only  limited  amount  of  information  on  the 
response  of  the  human  abdomen  to  impact  loads  from  steering 
wheels  [5,6].  Response  to  belt  loading  has  been  obtained 
essentially  from  testing  with  pigs  [2].  All  of  these  responses 
have  been  presented  in  the  form  of  force-deflection  curves  at 
various  impact  speeds.  The  force  is  measured  at  the  impactor  and 
the  deflection  is  the  total  deflection  of  the  abdominal  segment. 
The  DYNAMAN  model  was  set  up  so  that  the  necessary  disp  acemen 
and  forces  were  being  generated.  It  is  i not  P°^£le  ??e 
the  force  deflection  curve  directly  within  DYNAMAN,  and  the 
appropriate  dynamic  tables  were  imported  into  QuattroPro  to 
generate  the  necessary  curves  as  cross-plots. 
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OYNAMAN  VIEW 


TIME  CM5EC) 


Figure  3.  Model  of  a  belted  Hybrid  III  with  abdomen  component 


DYNAMAM  VIEW 
TIMECW5EC)  125. 


D.ao  14.73  ZM.45  MM.  IB  5B.S0 


Figure  4.  Model  of  Hybrid  III  with  abdomen  at  time  of  maximum 
penetration. 

A  number  of  different  approaches  were  taken  in  modelling  the 
spring-damper  system.  Instead  of  a  single  spring-damper,  a  set 
of  six  spring-dampers  were  defined  with  the  stiffness  equally 
distributed  amongst  them.  Such  a  model  was  built  to  verify  the 
overall  accuracy  of  the  original  model,  and  to  provide  a  more 
detailed  description  of  the  deformation  pattern  across  the  front 
of  the  abdomen.  Figure  5  shows  an  example  of  this  model. 
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Figure  5.  Multiple  spring  damper  model  at  time  of  maximum 
penetration 


Another  simulation  was  undertaken  by  modelling  the  lap  belt 
using  as  sequence  of  segments  connected  by  free  joints.  This 
allowed  the  analysis  of  the  response  due  to  the  extended  surface 
of  the  belt,  which  is  not  feasible  in  the  normal  harness  belt 
model.  Figure  6  shows  the  setup  of  such  a  belt  using  seven 
segments.  The  belt  surface  was  modeled  using  small  planes  which 
were  connected  by  free  pin  joints.  The  ends  of  the  surface  were 
connected  to  the  anchors  by  ordinary  harness  belt  segments. 


Figure  6.  Lap  belt  modelled  as  sequence  of  planes 


The  simulations  showed  that  the  results  when  the  belt  was 
modelled  using  planes  and  when  it  was  modelled  as  a  regular 
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harness  belt,  were  quite  comparable.  There  was  some  discrepancy 
between  the  model  using  a  single  spring-damper  and  that  using 
multiple  spring-dampers.  The  difference  appeared  to  arise  from 
the  multiple  contacts  which  were  generated  between  the  impacting 
surface  and  the  small  planes  that  were  defined. to  represent  the 
front  surfa.ce  of  the  bag  cover.  This  calculation  overestimated 
the  total  force  acting  on  the  impactor,  though  the  overall 
deflection  was  about  the  same  in  both  models.  It  was  felt  that 
the  model  using  a  single  spring-damper  would  be  accurate  enough 
for  deciding  on  the  range  of  foam  stiffness  that  would  be 
necessary  in  the  design.  Figure  7  is  a  comparison  of  the  results 
from  an  early  simulation  and  from  a  test  using  a  foam-filled  bag 
abdomen.  The  test  results  were  used  as  a  means  of  validating  the 
results  of  the  simulation,  before  extensive  design  changes  were 
performed.  As  part  of  the  analysis,  a  simple  analytic 
calculation  was  also  made  of  the  relative  contributions  of  the 
foam  and  the  airbag  to  the  stiffness  of  the  response.  The 
calculations  were  based  on  simple  energy  balance,  where  the 
initial  kinetic  energy  of  the  impactor  was  distributed  between 
the  deformation  energies  of  the  foam  and  the  bag.  The  foam 
deformation  energy  was  based  on  its  elastic  and  damping 
properties ,  while  the  airbag  deformation  energy  was  calculated 
using  the  first  law  of  thermodynamics  for  a  perfect  gas. 


Abdomen  Impact 

Impact  Velocity  =  15.05  ft/sec 


|  ■  Moonc/a  teat  ♦  08rr*jto<ton 


Figure  7.  Comparison  of  abdomen  simulation  and  test. 

The  different  kinds  of  design  changes  that  were  possible 
using  the  model  were: 

•  Volume  of  airbag 

•  Stiffness  of  foam 

•  Use  of  a  single  or  two  airbags 

•  Mass  of  cover 
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3.  Neck  Modelling 

A  new  neck  has  been  developed  by  the  Vehicle  Research  Test 
Center  (VRTC) ,  based  on  a  finite  element  model  developed  by  EASI, 
which  is  more  flexible  than  the  Hybrid  III  neck  [7].  The 
head/neck  system  has  been  tested  under  low  to  moderate 
acceleration  pulses  (8  -  15  G)  in  both  the  frontal  and  lateral 
directions  and  the  head  and  neck  angles  appear  to  come  close  to 
the  corridors  that  have  been  established  from  volunteer  testing 
T8  91.  Figure  8  shows  the  general  setup  of  the  new  neck.  It 
consists  of  five  aluminum  disks  separated  by  five  butyl  rubber 
disks.  The  rubber  disks  are  elliptical  in  shape  and  have 
variable  sizes  for  their  semiaxes. 


Figure  8.  Setup  of  new  neck. 

One  additional  feature  of  the  neck,  are  a  pair  of  cables 
that  are  attached  between  the  head  and  the  base  of  the  neck  and  a 
pair  of  spring  coils  that  are  attached  to  the  end  of  each  cable. 
DYNAMAN  was  used  to  develop  a  lumped  mass  model  of  the  aluminum 
disks  and  the  intermediate  rubber  disks.  The  neck  has  been 
developed  and  tested  as  a  separate  component,  attached  <?nly  to 
the  head.  As  part  of  an  effort  to  integrate  the  neck  with  the 
complete  dummy  and  properly  attach  the  neck  to  the  thoracic 
spine,  design  changes  which  included  changing  the  size  of  the 
disks  were  examined. 
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The  basic  model  of  the  neck  was  built  using  five  mass 
segments  to  represent  the  five  aluminum  disks.  The  aluminum 
disks  were,  for  practical  purposes,  rigid  segments  and  the 
bending  of  the  neck  was  generated  from  the  stiffness  of  the 
rubber  disks.  A  simple  pin  joint  was  defined  between  adjacent 
segments.  Two  different  models  were  developed  to  represent  the 
stiffness  of  the  neck  during  rotation.  The  first  model  was 
created  using  a  pair  of  spring-dampers  to  represent  each  rubber 
disk.  During  the  rotation  of  the  neck,  one  of  the  springs  would 
compress  and  one  would  extend  which  would  represent  the  actual 
compression  and  tension  on  the  two  sides  of  each  physical  disk. 

In  this  model,  the  joints  between  the  mass  segments  were  free 
segments  and  all  the  stiffness  was  carried  by  the  spring-dampers. 
The  second  model  was  based  on  assigning  appropriate  joint  torques 
at  the  joint,  without  the  recourse  to  spring-dampers. 

A  simple  calculation  procedure  was  used  to  come  up  with  a 
rough  estimate  of  the  properties  of  the  spring-dampers  from  the 
semiaxes  of  the  individual  disks  and  elastic  moduli  found  from 
compression  and  tension  tests  on  the  disks.  The  procedure 
assumed  a  distributed  compressive  stress  for  a  deformed  disk 
proportional  to  the  geometric  deformation  at  the  given  location 
and  then  integrating  over  the  area  of  the  disk,  to  come  up  with  a 
force-deflection  relation.  The  damping  coefficient  was  estimated 
based  on  the  responses  found  from  compression  tests  carried  out 
at  different  loading  rates.  Such  a  scheme  ignored  the 
contribution  of  the  shear  stress  and  this  was  modeled  by  assuming 
a  single  multiplicative  factor  which  could  be  used  to  fit  the 
response  at  one  acceleration  pulse  and  one  direction.  The 
comparison  of  responses  between  test  and  simulation,  for  other 
pulses  and  directions  would  provide  an  evaluation  of  the  efficacy 
of  the  model. 

The  joint  torque  function  used  in  the  joint  torque  model  was 
based  on  finite  element  modelling  that  had  been  carried  by  EASI, 
which  provided  the  joint  torques  as  a  function  of  angle  for  each 
of  the  disks.  It  was  found  that  the  joint  torque  model  was  not 
as  efficient  as  the  spring-damper  model  in  predicting  the 
head/neck  response  for  different  pulses  and  directions.  In 
addition  to  modelling  the  rubber  disks,  the  two  cables  were 
modelled  using  two  harness  belts.  One  end  of  each  belt  was 
anchored  to  the  head,  and  the  other  end  was  connected  to  a  small 
mass  segment,  which  in  turn  was  connected  by  a  spring  to  the  base 
of  the  neck.  Each  belt  passed  through  a  slip  ring  which  modeled 
an  eye  in  the  central  disk  through  which  the  cable  passed. 

Figure  9  shows  the  setup  of  the  neck  model  with  spring-dampers. 
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Figure  9.  Model  of  neck  using  spring- dampers 


As  mentioned  previously,  the  model  was  exercised  by  imposing 
frontal  pulses  of  maximum  magnitude  of  8  G  and  15  G  and  a  lateral 
pulse  of  7  G.  The  shapes  of  the  pulses  were  basically  half-sine 
waves.  The  behavior  of  the  model  under  extension  is  also  being 
examined,  though  the  human  response  data  under  such  loading  is 
not  as  thoroughly  documented.  The  output  that  was  examined,  were 
the  head  pitch  angle  and  the  neck  angle.  The  measurement  of  the 
neck  angle  poses  a  special  problem  in  these  tests,  since  the 
actual  neck  is  a  flexible  structure  with  variable  rotation  angles 
along  the  length  of  the  neck.  For  these  simulations,  a  standard 
location  relative  to  the  base  of  the  neck  was  selected  based  on 
the  length  of  the  rotation  axis  found  in  the  volunteer  tests  [8]. 
The  angle  made  by  a  vector  connecting  this  point  to  the  occipital 
condyle  (head/neck  joint)  was  taken  to  be  the  neck  angle.  The 
angle  was  computed  using  QuattroPro,  using  as  input  the  tables  o 
angular  displacement  of  the  selected  points  generated  by  DYNAMAN. 
Figure  10  shows  the  head/neck  system  at  the  point  of  maximum 
flexion.  Figure  11  shows  a  comparison  of  the  output  from  one  of 
these  simulations  with  the  volunteer  data  for  the  8  G  frontal 
pulse.  As  can  be  seen,  the  maximum  head  angle  come  fairly  close 
to  that  of  the  volunteer  average,  but  the  maximum  neck  angle  fall 
below. 

A  number  of  different  design  parameters  are  being  varied  in 
the  neck  model  in  order  to  understand  the  relative  sensitivity  of 
the  response  to  these  parameters  and  to  help  in  selecting  the 
parameters  that  can  be  successfully  integrated  into  a  revised 
neck  model.  The  basic  constraints  in  the  design  effort  are  the 
length  of  the  neck  which  should  be  similar  to  the  human  neck  (at 
least  for  the  portion  exposed  above  the  thorax)  and  the 
attachment  points  to  the  head  and  thorax. 
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Figure  10.  View  of  head/neck  model  at  maximum  flexion 


Figure  11. 


Comparison  of  head/neck  angles  from  volunteer 
tests  and  simulation 


The  following  areas  of  the  model  are  being  varied  to  examine 
their  influence  on  the  response: 

•  size  of  the  rubber  disks 

•  stiffness  of  the  rubber  disks 

•  number  of  rubber  disks 

•  stiffness  of  spring  attached  to  cable 

•  initial  preload  on  cable  spring 
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ATB-V 
GEB0D-IV.2 
VIEW-MI  &  IMAGE-1 


Louise  Obergefell 
Armstrong  Laboratory 
Wright-Patterson  Air  Force  Base 


ATB-V 


♦  Deformable  Element  Option 

♦  Water  Forces  Analysis  Capability 
(WAFAC) 

♦  Robotic  Joint  Actuators 

♦  Structured  Output  for  Graphics  Programs 

♦  36  and  15  msec  HIC 

4  Joint  Function  Use  with  Euler  Joints 
4  General  Code  Portability  Modifications 
4  New  Manuals 
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♦  Dr.  Hashem 
Ashrafiuon 

♦  Interconnected 
Deformable  & 
Rigid  Bodies 

♦  Uses  Natural 
Frequencies  & 
Normal  Modes 


Model  Element  Model 


Water  Forces  Analysis  Capability 


♦  Developed  for  Coast 
Guard  by  GESAC,  Inc. 

♦  Body  &  Personal 
Flotation  Device 
(PFD)  Response  to 
Water  Forces 

♦  Forces  Include 

-  Buoyancy 

-  Wave  Excitation 

-  Added-Mass 

-  Damping 

-  Drag  &  Lift 

♦  Still  Water  &  Wave 
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Robotic  Joint  Actuators 


♦  Active  Driving 
Feature 

♦  Joint  Torques 
Dependent  on 
Feedback 

♦  User-Defined 
Feedback 


2 


Structured  Output  for  Graphics 
Programs 


♦  Compatability  With  VIEW,  IMAGE,  & 
DYNAMAN  Postprocessor 

♦  Reduced  Size  of  ATB  Output  File 

4  Designed  to  Allow  Easy  Addition  of  New 
Graphics  Capabilities 
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ATB  Long-Term 


♦  New  Restraint  Belt  Model 

♦  DYNAMAN  Airbag 

-  Developed  by  GESAC 

+  Advanced  Contact  Capabilities 

-  Developed  by  GESAC 

♦  Harness  Belt  Improvements 


GEB0D-IV.2 


♦  Improved  Human  Joint  Properties 

♦  Advanced  Dynamic  Anthropomorphic 
Manikin  (ADAM)  Data  Set 

♦  Supplement  to  Manual 
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VIEW-1 1 1 


TIME  (MSEC)  300 


♦  Use  Structured  ATB 
Output 

♦  Final  Release 

-  Program  Will  Be 
Maintained  & 
Distributed 

-  All  Future  Graphic 
Improvements  Will  Be 
Made  To  IMAGE 

♦  New  Manuals 


IMAGE-1 

Interactively  Manipulated  ATB  Graphics 
Environment  (IMAGE) 


♦  Currently  Under 
Development 

♦  Resides  on  Silicon 
Graphics  Workstation 

♦  Uses  GL  Libraries 
(Open  GL  in  Future) 

♦  Solid  Model  Graphics 

♦  Easy  Graphical  Button 
&  Mouse  Interface 

♦  Users’  Manual 


Proceedings  of  the  1995  ATB  Model  Users’  Colloquium 


IMAGE 


♦  Zoom  &  Viewing  Angle  Controlled  Interactively 

♦  Play  Forwards,  Backwards,  or  Step  Frame 

♦  Colors,  Lighting,  &  Resolution  Are  Adjustable  via 
Control  File 

♦  Reads  Same  Output  File  as  VIEW 

♦  Saves  Images  in  rgb  Format 


Schedule 


♦  ATB-V  Fall  1996 

♦  GEBOD-IV.2  Winter  1995 

♦  VIEW-III  Spring  1996 

♦  IMAGE-I  l"7 
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GENERAL 
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DYNAMAN  IMPROVES  PRODUCTIVITY 


MODEL  SETUP 
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SETTING  UP  MODEL 
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CORRECT  THE  PULSE 
REFERENCES  NUMBERS 


ADDING  SEGMENTS 


LLIQ 
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OPTIONS  FILE  PICTURES  PLOTS  TABLES  EXIT 
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MODEL  SETUP 
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REPRODUCE  WELL  ON  VIDEO  WIT 
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•  SOFTWARE  MAINTENANCE  CHARGE 


LIST  OF  ATTENDEES 


Mr.  Prasenjit  Adhikari 
Renfroe  Engineering,  Inc. 

13045  West  Hwy  62 
Farmington,  AR  72730 
501/846-8000  Fax:  501/846-8002 

Dr.  Nabih  Alem 

U.S.  Army  Aeromedical  Research  Laboratory 
P.O.  Box  620577 
Fort  Rucker,  AL  36362-0577 
334/255-6892  Fax:  334/255-6937 

Mr.  Martin  Andries 

Air  Force  -  JPATS  Program  Office 

2121  Equestrian  Drive  #1B 

Miamisburg,  OH  45342 

513/436-3260 

Dr.  Hashem  Ashrafiuon 
Villanova  University 
Mechanical  Engineering 
800  Lancaster  Avenue 
Villanova,  PA  19085-1681 
610/519-7791  Fax:  610-519-7312 
hashrafiu@ucis .  vill .  edu 

Dr.  Xavier  Avula 
University  of  Missouri  -  Rolla 
231  Mechanical  Engineering  Bldg 
Rolla,  MO  65401 
314/341-4585  Fax:  314-364-3351 
avula@umr.edu 

Mr.  Lindley  Bark 
Simula  Government  Products 
10016  S.  51st  Street 
Phoenix,  AZ  85044 
602/730-4447  Fax:  602-893-8643 
lbark@enet.net 

Mr.  Scott  Barnes 
Ohio  State  University 
1087  Afton  Rd. 

Columbus,  OH  43221 
614/421-2803 

Dr.  Scott  Batterman 

Consultants  Associates,  Inc. 

1415  Route  70  East 
Suite  307 

Cherry  Hill,  NJ  08034 
609/795-3993  Fax:  609/795-4454 

Mr.  Paul  Bedewi 


George  Washington  University 
20322  Beechwood  Terr.  #303 
Ashbum,  VA  22011 
703/729-8367  Fax:  703-729-8359 

Ms.  Gina  Bertocci 
University  of  Pittsburgh  -  UPARC 
915  William  Pitt  Way 
Pittsburgh,  PA  15238 
412/826-3138  Fax:  412-826-3138 
ginaber@pitt.edu 

Mr.  David  Biss 

Automotive  Safety  Analysis 
15506  Avery  Ave. 

Suite  Number  One 
Rockville,  MD  20855 
301/294-7800  Fax:  301-294-0812 
71630. 2432@compuserve .  com 

Dr.  Bruce  M.  Bowman 
UMTRI  Biosciences 
2901  Baxter  Rd. 

Ann  Arbor,  MI  48103 
313/936-1106 

bmb@biosci  .umtri.umich .  edu 

Dr.  Frank  Buczek,  Jr. 

Advanced  Biomeics,  Inc. 

4927  Fayann  St. 

Orlando,  FL  32812 
407/384-7464Fax  :  407-384-7168 

Mr.  Arjaan  Buijk 
The  MacNeal-Schwendler  Corp 
26899  Northwestern  Highway 
Suite  400 

Southfield,  MI  48037 
810/208-3315  Fax:  810-208-0915 
arjaan.buijk@macsch.com 

Dr.  Russell  Burton 

AL/CF-CA 

2509  Kennedy  Circle 

Brooks  AFB,  TX  78235-5118 

210/536-2446  Fax:  210-536-2371 

russell .  burton@platinum  .brooks .  af .  mil 

Dr.  John  Cavanaugh 

Bioengineering  Center 
Wayne  State  University 
818  W.  Hancock 
Detroit,  MI  48202 
313/577-3916  Fax:  313-577-3916 
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Dr.  Philemon  (Phil)  Chan 
Jaycor 

9775  Towne  Centre  Dr. 

P.O.  Box  85154 

San  Diego,  CA  92186 

619/453-6580  Fax:  619-552-9172 

Mr.  Huaining  Cheng 

Systems  Research  Laboratories,  Inc. 

2800  Indian  Ripple  Rd. 

Dayton,  OH  45440 
513/255-3328 

hcheng@tweety.al.wpafb.af.mil 

Mr.  David  Chng 
Liability  Research 
55  Depot  Rd. 

Goleta,  CA  93117 
805/964-0676 

Mr.  Charles  Dickerson 
Arndt  &  Associates,  Ltd. 

2202  West  Huntington  Drive 
Tempe,  AZ  85282 
602/438-2004  Fax  :  602-438-0898 

Dr.  Kennerly  Digges 

FHWA/NHTSA 

National  Crash  Analysis  Center 

Virginia  Campus 

20101  Academic  Way 

Ashbum,  VA  22011-2604 

703/729-8363  Fax:  703/729-8359 

Mr.  Andrey  Dudkin 
Jaycor 

9775  Towne  Centre  Dr. 

P.O.  Box  85154 

San  Diego,  CA  92186 

619/453-6580  Fax:  619-552-9172 

Dr.  Janet  Dufek 

Department  of  Exercise  &  Movement  Science 

University  of  Oregon 

1240  University  of  Oregon 

Eugene,  OR  97403-1240 

503/346-3391  Fax:  503-346-2841 

j  dufek@oregon .  uoregon .  edu 

Mr.  Ed  Eveland 
AL/CFBV 
Building  441 
2610  Seventh  St. 

Wright-Patterson  AFB,  OH  45433-7901 
513/255-5963 

eeveland@tweety .  al .  wpafb .  af .  mil 

Mr.  Patrick  Fay 
Fay  Engineering  Corp 


5201  E.  48th  Ave. 

Denver  CO  80216 
303/333-5209  Fax:  303-329-0687 

Mr.  Richard  Fay 
Fay  Engineering  Corp. 

5201  E.  48th  Ave. 

Denver,  CO  80216 
303/333-5209  Fax:  303-329-0687 

Dr.  John  Fleck 
J&J  Technologies  Inc. 

92  Henning  Dr. 

Orchard  Park,  NY  14127 
716/662-4294 

Mr.  Steve  Forrest 
Liability  Research 
55  Depot  Rd. 

Goleta,  CA  93117 
805/964-0676  Fax:  805-964-7669 
lri@rain.org 

Mr.  Keith  Friedman 
Friedman  Research 
1513  Portesuello 
Santa  Barbara,  C A  93105 
805/569-7025 

Mr.  David  Furey 
Simula 

10016  S.  51st  Street 
Phoenix,  AZ  85044 
602/730-4445  Fax:  602-896-8643 

Mr.  Tom  Gardner 
Columbia  University 
Orthopedic  Laboratory 
630  W.  168th  St. 

New  York,  NY  10032 

Mr.  Wesley  Grimes 

Collision  Engineering  Associates,  Inc. 

P.O.  Box  31900 
Mesa,  AZ  85275-1900 
602/655-0399  Fax:  602-655-0693 
wgrhnes@host.yab.com 

Mr.  Charles  Griswold 
C.J.  Griswold,  Inc. 

5214  Clarendon  Crest  Court 
Bloomfield  Hills,  MI  48302-2607 
810/626-2470  Fax:  810-626-2470 

Dr.  Salvadore  Guccione 
Naval  Biodynamics  Laboratory,  Box  29407 
New  Orleans,  LA  70189-0407 
504/257-3962  Fax:  504-257-5456 
salg@tecnet.jcte.jsc.mil 
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Mr.  Herbert  M.  Guzman 

Biodynamic  Research  Corp. 

9901  IH  10  W.,  Suite  1000 
San  Antonio,  TX  78230 
210/691-0281  210-691-8823 
figaro.@adl.com 

Capt  Joel  Hagan 
AFIT/ENY 

Wright-Patterson  AFB  OH  45433 
DSN  785-6565  Ext  4401 
jhagan@afit.af.mil 

Ms.  Alena  Hagedorn 
Transportation  Research  Center,  Inc. 

P.O.  Box  B37 

East  Liberty,  OH  43319 

513/666-4511  Fax:  513-666-3590 

Mr.  Brian  Herbst 

Liability  Research 
55  Depot  Rd. 

Goleta,  CA  93117 
805/964-0676 

Ms.  Velvet  Hutson 
7335  Newell  Apt.  C 
Wichita,  KS  67212 
316/721-9490  Fax:  316-689-3175 
vhutson®  wsuhub .  uc .  twsu .  edu 

Mr.  Hirotoshi  Ishikawa 

Japan  Automobile  Research  Institute  Inc. 

2530  Karima  Tsukuba  City 

Ibaraki  305 

Japan 

+  81  298  56  1111  Fax:  +81  298  56  1135 

Dr.  Ints  Kaleps 
AL/CFBV 
Building  441 
2610  Seventh  St. 

Wright-Patterson  AFB,  OH  45433-7901 
513/255-3665 

ikaleps@tweety .  al .  wpafb .  af .  mil 
Mr.  Alva  Karl 

Systems  Research  Laboratories,  Inc 
2800  Indian  Ripple  Rd. 

Dayton,  OH  45440 
513/255-3328 

Mr.  Scott  Kebschull 
Dynamic  Research,  Inc. 

355  Van  Ness  Ave.  #200 
Torrance,  CA  90501 
310/212-5211  Fax:  310-212-5046 
keb  1  @  ix .  netcom .  com 


Dr.  Tara  Khatua 

Failure  Analysis  Associates 
149  Commonwealth  Drive 
Menlo  Park  CA  94025 
415/326-9400  Fax:  415-328-2981 
kha@fail.com 

Mr.  Kerry  Knapp 

3632  N.  Schevene 
Flagstaff,  AZ  86004 
520/523-6771  Fax:  520-779-9621 

Mr.  Richard  Lawrence 
SEA,  Inc. 

7349  Worthington-Galena  Rd. 
Columbus,  OH  43085 
614/888-4160 

Ms.  Laura  Liptai 

University  of  California  -  Davis 
Biomedical  Engineering 
514  Central  Avenue 
Alameda,  CA  94501 
510-376-1240  Fax:  510-376-1245 
liptai®  venus .  engr .  ucdavis  .edu 

Dr.  Deren  Ma 

Systems  Research  Laboratories,  Inc 
2800  Indian  Ripple  Rd.,  WP441 
Dayton,  OH  45440 
513/255-3328 

dma@tweety .  al .  wpafb .  af .  mil 

Mr.  Jeffrey  Marcus 

FAA-CAMI 

AAM-630 

P.O.  Box  25082 

Oklahoma  City,  OK  73125 

405/954-5555 

Mr.  Joseph  McCarthy 
Becker  Transportation  Safety  Inc. 
1725  North  Talbot  Rd.  R.R.  #1 
Windsor,  Ontario  N9A  6J3 
519/737-7233  Fax:  519-737-7203 

Mr.  Joe  McEntire 

USAARL 

MCMR-UAD-CI 

Fort  Rucker,  AL  36362-0577 

334/255-689 6  Fax:334/255-6937 

Mr.  Brian  McHenry 
McHenry  Consultants,  Inc. 

103  Brady  Court,  Ste  200 
CaryNC  27511 
919/469-3310 
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Mr.  Peter  Mendoza 
MacNeal-Schwendler  Corp. 

815  Colorado  Blvd. 

Los  Angeles,  CA  90041 
213/259-3860  Fax:  213/2594992 
peter .  mendoza@macsch.com 

Dr.  Thomas  Moore 
USAF  -  AL/CFB 
Building  441 
2610  Seventh  St 

Wright-Patterson  AFB,  OH  45433-7901 
513/255-3603 

Ms.  Ellen  Moore 

Systems  Research  Laboratories,  Inc 

2800  Indian  Ripple  Rd. 

Dayton,  OH  45440 
513/255-3328 

Mr.  Lam  Nguyen 
R.D.  Jablonsky  Inc. 

2059  E.  Foothill  Blvd. 

Pasadena,  CA  91107 
818/405-8997 

Mr.  Gert  Nilson 

GM  R&D  Center 

Auto  Safety  &  Health  Research 

30500  Mound  Road 

Warren,  MI  48090-9055 

810/696-5218  Fax:  810-696-5150 

ifggen@troy .  ifa .  gmeds .  com 

Mr.  Ronald  Nordhagen 
Collision  Safety  Engineering,  Inc. 

150  South  Moutainway  Drive 
Orem,  UT  84058 
801/229-6252 

Dr.  Louise  Obergefell 
AL/CFBV 
Building  441 
2610  Seventh  St. 

Wright-Patterson  AFB,  OH  45433-7901 
513/255-3665 

lobergef@tweety .  al .  wpafb .  af .  mil 
Dr.  Peter  Orner 

Internal  Medicine  &  Biomechanics  of  Injury 
12910  Lomas  Verdes  Dr. 

Poway,  CA  92064-1250 
619/487-1855  Fax:  619487-9743 

Mr.  Joe  Pellettiere 
AL/CFBV 
Building  441 
2610  Seventh  St. 

Wright-Patterson  AFB,  OH  45433-7901 


513/255-5963 

j  pellet@tweety .  al .  wpafb .  af .  mil 

Mr.  Chris  Perry 
USAF  -  AL/CFBE 
Building  824 

Wright-Patterson  AFB,  OH  45433 
513/255-3122 

Mr.  Kenneth  Potempa 
Air  Force,  HSC/XRS 
2510  Kennedy  Circle,  STE  220 
Brooks  AFB,  TX  78235-5120 
210/5364452 

Mr.  Howard  B.  Pritz 
2416  Fishinger  Rd. 

Columbus,  OH  43221 
614/451-6898 

Mr.  John  Quartuccio 
Naval  Air  Warfare  Center 
P.O.  Box  5152 
Code  4.6. C  MS-15 
Warminster,  PA  18974-0591 
215/441-7602  Fax:  215/441-3765 

CAPT  Edward  Rivers 
USAF  -  AL/SDN 

Wright-Patterson  AFB,  OH  45433 
513/255-7619 

Ms.  Annette  Rizer 

Systems  Research  Laboratories,  Inc. 

2800  Indian  Ripple  Rd. 

Dayton,  OH  45440 
513/255-3328 

arizer@tweety.al.wpafb.af.mil 

Ms.  Marjorie  Seeman 
Naval  Biodynamics  Laboratory 
Box  29407 

New  Orleans,  LA  70189-0407 
504/257-3962  Fax:  504-257-5456 

Mr.  Brian  Self 
AL/CFBV 
4308  S.  Abby  Ct. 

Salt  Lake  City,  UT  84123 
801/2694016 

Dr.  Tariq  Shams 
GESAC,  Inc. 

Route  2  Box  339A 
Keameysville,  WV  25430 
304/267-6849  Fax:  304-267-6821 
gesac@access .  mountain .  net 
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Mr.  Michael  Shanahan 

Becker  Transportation  Safety  Inc. 

1725  North  Talbot  Rd.  R.R.  #1 
Windsor,  Ontario  N9A  6J3 
519/737-7233  Fax:  519-737-7203 

Dr.  Ed  Sieveka 

University  of  Virginia 

Dept,  of  Mechanical  Engineering 

Thornton  Hall 

Charlottesville,  VA  22903 

804/924-7075  Fax:  804-982-2037 

ems@virginia.edu 

Ms.  Jeanne  Smith 

Systems  Research  Laboratories,  Inc. 

2800  Indian  Ripple  Rd. 

Dayton,  OH  45440 
513/255-3328 

j  smith@tweety .  al .  wpafb .  af .  mil 

Mr.  Manohar  Srinivasan 
General  Motors  Corp. 

7th  Floor  New  Center  One 
3031  W.  Grand  Blvd. 

Detroit  MI  48232 
313/974-1369Fax  :  313-974-0690 

Dr.  Chi-Ming  Tang 

State  University  of  New  York  -  Geneseo 
Department  of  Mathematics 
1  College  Circle 
Geneseo,  NY  14454 
716/245-5386  Fax:  716-245-5005 
tang@uno .  cc .  geneseo .  edu 

Mr.  David  Tung 
Delphi  Automotive 
8801  Shadycreek  Drive 
Centerville,  OH  45458 
513/356-2338  Fax:  513-356-2303 

Mr.  Ruben  van  Schalkwijk 
TNO  Crash-Safety  Research  Centre 
Schoemakerstraat  97 
2628  UK  Delft 
Holland 

+31  15  696951  Fax:  +31  1562  4321 
vanschalkwij  k@  wt  .tno.nl 

Mr.  Michael  Varat 
JFK  Engineering 
5636  La  Cumbre  Rd. 

Somis,  CA  93066 
805/386-3388  Fax:  805-386-4892 

Mr.  Liming  Voo 

Medical  College  of  Wisconsin 


Research  151 

5000  W.  National 

Milwaukee,  WI  53295 

414/384-2000  Ext.  1529  Fax:  414-382-5374 

1  voo@post .  its .  mew .  edu 

Mr.  Don  Weerappuli 

Ford  Motor  Co. 

20000  Rotunda  POB  2053 
Mail  Drop  66 
Dearborn,  MI  48121-2053 
313/248-9953  Fax:  313-248-7847 

Mr.  Greg  Weisenfeld 
Honda  R  &  D 
1900  Harpers  Way 
Torrance  CA  90501 
310/781-5844 

Dr.  Mariusz  Ziejewski 
North  Dakota  State  University 
2363  20th  Avenue  South 
Fargo,  ND  58103 
701/232-9223  Fax:  701-232-9223 
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